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ABSTRACT 

Using high spatial resolution Hubble Space Telescope Wide Field Camera 3 and Advance Camera for 
Surverys imaging of resolved stellar populations, we constrain the contribution of thermally-pulsing 
asymptotic giant branch (TP-AGB) stars and red helium burning (RHeB) stars to the 1.6 /im near- 
infrared (NIR) luminosities of 23 nearby galaxies, including dwarfs and spirals. The TP-AGB phase 
contributes as much as 17% of the integrated F160VK flux, even when the red giant branch is well 
populated. The RHeB population contribution can match or even exceed the TP-AGB contribution, 
providing as much as 21% (18% after a statistical correction for foreground) of the integrated F160iy 
light. We estimate that these two short lived phases may account for up to 70% of the rest-frame 
NIR flux at higher redshift. The NIR mass-to- light (M/L) ratio should therefore be expected to vary 
significantly due to fiuctuations in the star formation rate over timcscalcs from 25 Myr to several 
Gyr, an effect that may be responsible for some of the lingering scatter in NIR galaxy scaling relations 
such as the TuUy-Fisher and metallicity-luminosity relations. We compare our observational results to 
predictions based on optically derived star formation histories and stellar population synthesis (SPS) 
models, including models based on the 2008 Padova isochrones (used in popular SPS programs) and 
the updated 2010 Padova isochrones, which shorten the lifetimes of low-mass (old) low-metallicity 
TP-AGB populations. The updated (2010) SPS models generally reproduce the expected numbers of 
TP-AGB stars in the sample; indeed, for 65% of the galaxies, the discrepancy between modeled and 
observed numbers is smaller than the measurement uncertainties. The weighted mean model/data 
number ratio for TP-AGB stars is 1.5 (1.4 with outliers removed) with a standard deviation of 0.5. 
The same SPS models, however, give a larger discrepancy in the _F160M^ flux contribution from the 
TP-AGB stars, over-predicting the flux by a weighted mean factor of 2.3 (2.2 with outliers removed) 
with a standard deviation of 0.8. This larger offset is driven by the prediction of modest numbers of 
high luminosity TP-AGB stars at young (< 300 Myrs) ages. The best-fit SPS models simultaneously 
tend to under-predict the numbers and fluxes of stars on the RHeB sequence, typically by a factor of 
2.0 ± 0.6 for galaxies with significant numbers of RHeBs. Possible explanations for both the TP-AGB 
and RHeB model results include: (1) difficulties with measuring the SFHs of galaxies especially on the 
short timescales over which these stars evolve (several Myrs); (2) issues with the way the SPS codes 
populate the CMDs (e.g. how they handle pulsations or self extinction), and/or (3) lingering issues 
with the lifetimes of these stars in the stellar evolution codes. Coincidentally these two competing 
discrepancies — over-prediction of the TP-AGB and under-prediction of the RHeBs — result in a 
predicted NIR M/L ratio largely unchanged for a rapid star formation rate, after correcting for these 
effects. However, the NIR-to-optical flux ratio of galaxies could be significantly smaller than AGB-rich 
models would predict, an outcome that has been observed in some intermediate redshift post-starburst 
galaxies. 

Subject headings: galaxies: stellar content — stars: AGB and post-AGB — stars: Hertzsprung-Russell 
diagram — galaxies: fundamental parameters 
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1. INTRODUCTION 



One of the primary objectives of extragalactic obser- 
vational astronomy is to measure and track the growth 
of stellar m ass in ga laxies across cosmic time (e.g . 
Bundv et al.l[20 05: Font ana et al.ll2006l:IIIbCTt et al.ll201C)l: 



Pozzetti et alBoiQ : ,Vulcani et al.H2010D . To accomplish 



this task, rest-frame ultra-violet (UV) through near- 
infrared (NIR) observations have been obtained for hun- 
dreds of thousand s of galaxies (e.g.fGiav alisco et al. 2004; 
iDavis et al.ll2007t iSanders et al.ll2007[r However, the in- 
terpretation of these observations requires stellar popula- 
tion synthesis codes that incorporate models of the initial 
mass function, star formation histories, and stellar evolu- 
tion tracks (e.g. Bruzual & Chariot 2003). The detailed 
prescriptions for these inputs can affect the resulting es- 
timates of stellar population age, and to tal stellar mass 
(jMaraston et al.|[200a ITlbert et al.llMol) . 

Until recently, NIR passbands were assumed to provide 
an ideal window o n the stellar masses of galaxies (e.g. 
iBundv et al.l l200"5'). Compared to optical and UV pass- 
bands, NIR wavelengths are significantly less affected by 
massive main sequence stars formed in bursts of star for- 
mation, which can decrease the mass-to- light (M/L) ratio 
in the optical passbands. NIR observations are also less 
affected by dust obscuration, which can increase the M/L 
ratio at shorter wavelengths. In addition, deep Spitzer 
IRAC observations provide an ideal window on the rest- 
frame NIR fluxes of high redshift galaxies, and have been 
used extensively for estimating stellar masses. 

Unfortunately, while massive main sequence stars do 
not have a large impact on the NIR luminosities of galax- 
ies, intermediate- mass (2 — 10 Mq) evolved stars have 
been shown to contribute significantly to integrated NIR 
fluxes, even when they represent a negligible contribu- 
tion t o the stellar mass (jPersson et al.|[l983l : iFrogel et al.l 
|1990[) . In recent years, renewed effort has been given 
to understanding the contribution of thermally-pulsing 
asymptotic giant branch (TR^AGB^ s tars t o the NIR 
M/L ratios of galaxies ( M araston et all 120061) . The TP- 
AGB represents a brief period (1 — 2 Myr) of double 
shell burning at the end of stellar evolution. During this 
phase, a star swells, undergoes pulsations, and ultimately 
loses as much a s 80% of its stellar mass before fading to 
a white dwarf ( Ibcn fc Rcnz ini|[l983t iVassiliad is fc WoodI 
[1993; Kcnnicutt-el_aL _1994|). The most massive of these 
TP-AGB stars can be very luminous in the NIR, ex- 
ceeding the luminosity of the tip of the red giant branch 
(TRGB) by several magnitudes. Models that neglect TP- 
AGB stars have been shown to over-estimate the masses 
of distant galaxies by factors of t wo or more in com pari- 
son to models that include them (jllbert et al.ll2010[ ). 

TP-AGB stars are now routinely included in popu- 
lation synthesis models of galaxies, although in differ- 
ent proportions depending on the technique adopted 
(see Chariot & Bruzual 1991; Brcssan et al. 199; 



Maraston et al.l 120061 : iBruzual Al 120071 : IConrov et al 
2009). Many issues still remain, primarily because late 
stage stellar evolution is difficult to follow from first prin- 
ciples. Stellar evolution codes require knowledge of hard 
to model processes, such as: recurrent third dredge-up 
events, hot-bottom burning, long period variability, and 
mass loss (Marigo & Girardi 2007). To account for these 
processes, modelers often resort to simplified TP-AGB 



stellar evolution codes with parame ters tuned to observa- 
tional data sets (jMarigo et al.ll2008l) . Currently, the most 
complete data sets of evolved stars come from studies 
of the Large and S mall Magellanic C louds (iFroeel et al 
T990t iCioni et al.l [1991 iBlum et al.l [2001 lBover"era! 



20m . The TP-AGB evolution code used to build the 
Padova isochrones (Girardi & Marigo 2007; Marigo et 
al. 2008) has been successfully tuned to reproduce the 
numbers and optical luminosities of TP-AGB stars in 
these systems. 

Unfortunately, the stellar populations of the Magel- 
lanic Clouds only span a narrow region of age and metal- 
licity. Codes tuned to the Magellanic clouds can fail 
dramatically when used to predict the resolved stel- 
lar popula tions o f other nearby galaxies. For instance, 
Gullicuszik et al.l ([2008b) showed that these codes over- 
predicted the numbers of carbon-rich AGB stars in the 
Leo II dSph by a factor of six. Similarly, in nearby 
(2.5 Mpc) dwarf irregular galaxy KKH 98, the models 
have been found to over-predict the numbers of TP- 
AGB stars by factors of 2-3 compared to observations 
(i Melbourne et al.ll2010( ). These codes also had difficulty 
modelling more massive metal-rich galaxies at larger dis- 
tances. For instance, the Virgo Cluster shows a de- 
ficiency of AGB star s compared to model predictions 
([Williams et al.|[2007l) , and a sample of intermediate red- 
shift post-starburst galaxies shows spectral energy distri- 
butions tha t rule out large fl ux contributions from TP- 
AGB stars (jKriek et al.ll2010(l . 

Star count analysis has been extended to galaxies be- 
yond the Local Group, where resolved stellar popula- 
tions are best studied with HST (e.g. the ACS Ne arby 
Galaxy Survey Treasury, ANGST lDalcanton et aLl2 009'). 
or with adaptive optics on large ground based tcle- 
scopes ( GuUieuszik et al. 2008a; Melbourne et al. 201(3; 



Davidgell201(|) . In one of the largest such studies to date, 
Girardi et aLl ([2010D found that that the 2008 Padova 



isochrones were over-predicting the TP-AGB in optical 
HST observations of 10 old, metal-poor galaxies from 
the ANGST sample. However, the models could be 
brought into agreement with the data by lowering the 
estimated lifetimes of low- mass (old) , low-metallicity TP- 
AGB stars, making them roughly equivalent to the life- 
times of higher mass (younger) TP-AGB stars. These 
revisions have been incorporated into the 2010 versions 
of Padova stellar evolution codes and isochrones [IB- 
While significant effort has recently been given to the 
TP-AGB phase, another phase may be equally impor- 
tant for understanding NIR M/L ratios of galaxies. Lu- 
minous red helium burning stars (RHeBs) are massive 
(i.e. > 3.5Mq) core helium burning stars that form a 
tight se quence at the luminous end of optical and NIR 
CMD s (|Dohm-Palmer fc Skillm"m] 120021 : iMcQuinn et~all 
120 lit ). These stars have even shorter lifetimes than TP- 
AGB stars, and the RHeB sequence is only obvious in 
CMDs with a burst of star formation younger than ~ 300 
Myrs. As with TP-AGB stars, these stars are difficult 
to model, with significant uncertainties associated with 
convective/mixing processes (overshooting and rotation) 
and mass loss. Unfortunately, less attention has been 
given to this phase of stellar evolution for stellar popula- 
tion synthesis models, even though at high redshift where 

http:/stcv. oapd.inaf.it/cmd 
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Figure 1. Example CMDs of the stars in program galaxy UGC-4305-1. These CMDs contain only stars that are matched across all four 
optical and NIR bands; hence, stars that are blue and faint tend to be missing from the optical CMD (upper left). Note, the optical and 
optical/NIR hybrid (upper right) CMDs provide a much larger color spread for IR luminous stars, compared with the NIR only CMD 
(bottom left). This color spread makes it easier to distinguish between RHcB stars and TP-AGB stars. We therefore use the optical/NIR 
hybrid CMD to identify sequences of evolved stars (see Figure [2J. 



there is significant ongoing star formation, RHeBs will be 
one of the dominant contr ibutors to the NIR luminosity 
(see iDalcanton et"aLll2012D . 

In this paper, we build on IGirardi etall (poTol ). now 
examining the TP-AGB and RHeB stars within a di- 
verse sample of 23 dwarf and spiral galaxies in the 
nearby universe, many with significant on-going star for- 
mation. Individual stars within these galaxies are re- 
solved with high spatial resolution Hubble Space Tele- 
scope (HST) observations in the optical (ACS) and NIR 
(WFC3). HST crowded- field photometry techniques 
provide the distributions of stars in color-magnitude 
space. The ACS observations are deep enough to con- 
strai n the star formation his t ories ( SFHs) of the s e galax - 
ies (IWilliams et all I2009bl IMot IWeisz et all [20lTal) . 
while the WFC3 data provide constraints on the NIR lu- 
minosities of the TP-AGB, RHeB, a nd red giant branch 
(RGB) stars ([Dalcanton et al.ll2012f ). First we calculate 
the fraction of the 1.6 /im galaxy flux contributed by 
TP-AGB and RHeB stars as a function of the population 
age, a number that can be used to correct the NIR M/L 
ratios of galaxies for evolved stellar populations. Then 
we compare the observations to the numbers and fluxes 
of TP-AGB and RHeB stars predicted by stellar pop- 
ulation synthesis (SPS) models and the 2008 and 2010 
Padova isochrones. Finally we discuss the implications 
of IR luminous stars for models of high redshift galax- 



ies. The next paper in this series Rosenfield et al. (in 
preparation) will explore these results further and pro- 
vide updates to the Padova stellar evolution codes where 
needed. 

This paper is arranged as follows. Section 2 describes 
the optical and NIR HST observations of the sample 
galaxies. Section 3 provides the flux fractions con- 
tributed by TP-AGB and RHeB stars from both our 
data and simulations. Section 4 examines the results 
in more detail and discusses the implications of studies 
at high redshift. Section 5 summarizes our conclusions. 
Magnitudes are reported in the Vega system, and we as- 
sume the can onical ACDM cosm ology with f2j\/ = 0.3 
and = 0.7 HSpergel et all [20071 ). 

2. THE DATA: RESOLVED STELLAR POPULATIONS 
FROM HST 

Studies of resolved stellar populations require very high 
spatial resolution (< 0.1") imaging for galaxies outside 
the Local Group . Even at the resolution of HST, these 
studies are only possible within roughly the local 4 Mpc 
volume. ANGST provided the flrst uniform observa- 
tional data-set of optical HST imaging of galaxies within 
the local volume. A subset of the ANGST sample has 
subsequently been observed wit h high spatial resolutio n 
NIR imaging with HST WFC3 (iDalcanton et aLllMl) . 
In this paper, we use the optical and NIR HST ANGST 
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Table 1 

Properties of the Sample Galaxies Measured from the SFH Fitting Routine CalcSFH 



Galaxy Distance Metallicity Fraction of Mass Fraction of Mass 

Modulus [M/H] Younger than 2 Gyr Younger than 0.3 Gyrs 



DD071 


27 


67 


-1 


09 


± 





31 


1 


43C-02 


± 


6 


Oe 


03 





OOe-00 


± 





Oe-00 


DD078 


27 


84 


-1 


15 


± 





16 


1 


88e-02 


± 


3 


9e 


03 


1 


OOe-04 


± 


1 


Oe-04 


DD082 


28 


04 


-1 


11 


± 





16 


1 


75e-02 


± 


3 


6e 


03 


7 


OOe-04 


± 


1 


Oe-04 


ESO540-030 


27 


76 


-1 


06 


± 





24 


6 


15e-02 


± 


2 


4e 


02 


1 


30e-03 


± 


6 


Oe-04 


HS117 


27 


93 


-0 


61 


± 





30 


5 


30e-02 


± 


5 


5e 


03 


1 


60e-03 


± 


1 


Oe-04 


IC2574-SGS 


27 


98 


-0 


97 


± 





43 


9 


27e-02 


± 


8 


7e 


03 


1 


44e-02 


± 


5 


le-03 


KDG73 


27 


90 


-1 


28 


± 





14 


6 


36e-02 


± 


1 


5e 


02 


9 


70e-03 


± 


4 


5e-03 


KKH37 


27 


66 


-1 


04 


± 





11 


2 


55e-02 


± 


3 


7e 


03 


2 


lOe-03 


± 


6 


Oe-04 


M81-DEEP 


27 


78 


-0 


41 


± 





28 


2 


09e-02 


± 


1 


le 


02 


2 


90e-03 


± 


1 


8e-03 


NGC0300-WIDE1 


26 


55 


-0 


75 


± 





36 


8 


83e-02 


± 


2 


Oe 


02 


1 


96e-02 


± 


6 


5e-03 


NGC2403-HALO-6 


27 


51 


-0 


65 


± 





51 


9 


88e-02 


± 


6 


5e 


03 


1 


17e-02 


± 


4 


le-03 


NGC2976-DEEP 


27 


73 


-0 


78 


± 





40 


1 


80e-03 


± 


4 


7e 


03 


2 


OOe-04 


± 


1 


Oe-04 


NGC3077-PHOENIX 


27 


95 


-0 


94 


± 





11 


5 


lOe-03 


± 


6 


Oe 


04 


3 


80e-03 


± 


6 


Oe-04 


NGC3741 


27 


49 


-1 


27 


± 





16 


1 


lle-01 


± 


1 


le 


02 


2 


06e-02 


± 


4 


6e-03 


NGC4163 


27 


36 
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19 


± 
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± 


4 


3e 


03 


3 


50e-03 


± 


7 


Oe-04 


NGC7793-HALO-6 


27 


91 


-0 


69 


± 





34 


6 


07e-02 


± 


9 


le 


03 


1 


19e-02 


± 


1 


8e-03 


SCL-DEl 


28 


22 


-1 


19 


± 





17 


2 


36e-02 


± 


1 


le 


02 


4 


OOe-04 


± 


3 


Oe-04 


UGC4305-1 


27 


64 


-1 


15 


± 





12 


1 


03e-01 


± 


9 


2e 


03 


2 


39e-02 


± 


5 


le-03 


UGC4305-2 


27 


64 


-1 


14 


± 





19 


1 


lle-01 


± 


1 


le 


02 


2 


58e-02 


± 


6 


4e-03 


UGC4459 


27 


79 


-1 


16 


± 





16 


6 


60e-02 


± 


2 


le 


02 


9 


OOe-03 


± 


1 


3e-03 


UGC5139 


27 


91 


-0 


78 


± 





36 


1 


03e-01 


± 


1 


4e 


02 


1 


97e-02 


± 


5 


6e-03 


UGC8508 


27 


04 


-1 


25 


± 





14 


5 


62e-02 


± 


1 


2e 


02 


1 


33e-02 


± 


2 


9e-03 


UGCA292 


27 


54 


-1 


51 


± 





06 


1 


97e-01 


± 


3 


5e 


02 


4 


96e-02 


± 


1 


7e-02 



observations of 23 nearby galaxies to study luminous TP- 
AGB and RHeB stars. 

2.1. Optical ACS Imaging and Star Formation 
Histories 

Deep multi-band optical HST imaging of a complete 
set of galaxies within the local 4 Mpc volume was either 
obtained by the ANGST team or culled from the archive 
and included in the ANGST program. ANGST galaxies 
were observed in at least two filters to provide color and 
luminosity information for the stars. Each galaxy was 
observed in the red F814W filter. At least one bluer 
band was also obtained, usually in the F475W, i^555W^, 
or F606W bands. 

The basic image reductions were described in 
iDalcanton et al.l ()2009( ). Photometry of these fields were 
obtained with the HST crowded-field ph otometry pack - 
age DOLPHOT, a version of HSTPHOT (|Dolphinll2000l ). 
which has been optimized for use with ACS and WFC3. 
The sensitivit ies of thes e photometric data-sets are pro- 
vided in Dalca nton et al . (2009) , but typically were deep 
enough to reach the red clump and the main sequence 
turn-off for populations younger than 1 Gyr. Figure [T] 
shows an example optical CMD. 

The optical multi-band photometry was used to con- 
strain the SFH of each galaxy. The numbers and po- 
sitions of stars across color-magnitude space are set by 
stellar evolution and the SFH of each galaxy. Both the 
youngest and oldest stellar populations are thought to 
be well const rained by the CMD s in the ANGST sample 
(see tests in iWeisz et al.l [2011 al) . The youngest popu- 
lations are constrained by luminous main sequence stars 
and evolving supergiant stars. Older populations are well 
constrained by the RGB, which becomes well-populated 
for galaxies older than ~ 2 Gyrs. 

The globa l SFHs of the ANG ST dwarf galaxies are 
described in IWeisz et al.l (|2011aD , and were determined 
with the 2-D CMD fitting routine CalcSFH (part of the 



MATCH package jDolEhi^ [2000) . We constructed the 
SFH of each galaxy based on the numbers of stars within 
color magnitude bins on the observed CMD, with color 
bins of size 0.05 mags, and magnitude bins of 0.1 mags. 
We used 71 logarithmic time bins from 4 Myrs to 12 
Gyrs old, and 24 different metallicity bins ranging from 
metallicities of [M/H] — —2.25 to 0.05. We assumed 
a single-slope power-law IMF with a spectral index of 
-1.30 over a mass range of 0.1 — 120 Mq, and a binary 
fraction of 0.35 with a flat secondary mass distribution. 
The difference between our selected IMF and a Kroupa 
IMF (Kroupa 2001) is negligible, as the ANGST CMDs 
are limited to stellar masses > 0.8 Mq. 

The methods for es timating uncertaint ies in the SFHs 
are also described in IWeisz et"aLl (j2011aD . We used the 
standard MATCH routine Monte Carlo approach to es- 
timate both the random and systematic uncertainties. 
For each Monte Carlo run, the observed CMD was ran- 
domly resampled and refit with CalcSFH. Additive errors 
in Mboi and log(Te//) were introduced when generating 
the model CMDs for these solutions. This method was 
developed to account for the full range of systematic dif- 
ferences between isochrone sets that use different pre- 
scriptions for various phenomena such as stellar rotation 
and convective overshooting. 

For this paper, we focus on the sub-regions of the 
ANGST galaxies that were observed in both the optical 
with ACS and the NIR with WFC3. These sub-regions 
are smaller than the full ACS fields, and may have dif- 
ferent SFHs than the full fields. We therefore re-ran 
the SFH modehng codes using the optical photometry 
of just the overlap region, but following the procedure 
and binning used in IWeisz et al.l (j2011a ). This reanaly- 
sis typically did not change the relative amount of star 
formation in each time bin but only the overall scaling. 
Table[T]gives the distance modulus, metallicity, and frac- 
tions of young stars for each galaxy as measured by the 
SFH routine CalcSFH. Table [T] also provides uncertain- 
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Figure 2. CMDs of the program galaxies in tlie (_F814VF — -F16014^) bands (top left). Regions that identify different features in the 
CMDs are shown: RHeB (blue), TP-AGB (red), and RGB (brown). Model CMDs created from the best- fit star formation histories are 
also shown for the SPS codes based on the 2008 Padova isochrones (top middle) and 2010 isochrones (top right). The models do a good 
job of reproducing the RGB. As expected the 2008 models significantly over-predict the numbers of TP-AGB stars (red points), while the 
2010 models are a much better match. However, both models tend to over-predict the luminous TP-AGB populations. The SPS models 
also tend to under-predict the numbers of RHeB stars (blue points), especially at the luminous end. 

respectively. The observations and image reduction o f 
the WFC3 data are described in lDalcanton et all (|2012l ). 
As with the optical data, photometry of the WFC3 ob- 
servations was done with the DOLPHOT package, which 
has been updated to include a module for the processing 
of WFC3 data. 

Figure [T] shows an example NIR CMD produced from 
these data. Typical uncertainties for the stellar photom- 
etry range from 0.01 mags at the FY'oQW ~ 18 to 0.10 
mags at i^lGOVF ~ 24. Each galaxy is observed to sev- 
eral magnitudes below the TRGB. These limits are not 
typically faint enough to detect the main sequence turn 
off for stars older than a couple hundred Myrs, or the 
well-populated red clump. However, the CMDs do show 
AGB, RGB, and RHeB populations (Figure 1). 



ties for these parameters. Both metallicity and fractions 
of young stars are characterized by relatively large un- 
certainties compared to Local Group galaxies for which 
the SFH has been determined with the same method 
jDolphin et al."2005'; 'Hohzman et al." 2006; G allart et all 
[foOS; WiUiams et al. 2009a; Wcisz ct al. 201113). As dis- 
cussed in IGirardi et al.l (j201Ql , larger SFH uncertainties 
are the price to pay for observing large samples of TP- 
AGB stars with only a single RST pointing per galaxy, 
and with minimal contamination from foreground stars. 

2.2. HST NIR WFC3 Imaging 

HST NIR images of a subset of the ANGST sample 
were obtained during Cycle 17, in program SNAP-11719. 
Imaging was obtained in both the FllOM^ and Fl&OW 
filters, with total exposure times of 597.7 s and 897.7 s 
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Figure 2. continued 

2.3. Catalogue Matching 

By design there is significant overlap between the opti- 
cal and NIR images of each galaxy. We generate optical 
through NIR matched catalogues to identify TP-AGB 
and RHeB stars. To do a proper transformation between 
the two coordinate systems, we select ~ 150 stars that 
are bright in both the optical and NIR data sets and 
that spatially span the entire overlap region between the 
WFC3 and ACS images. Starting with the optical and 
NIR catalogues from DOLPHOT, we cull the lists to only 
include stars that are in the spatial overlap region. Then 
we select all of the reasonably bright red stars in each 
dataset, with optical color > 0.7 mags and F814W^ < 26 
mags, and IR color (FllOW^ - FimW) > 0.5 mags and 
F160W < 24 mags. We sort these two lists by luminos- 
ity and select 150 stars roughly evenly spaced across the 
image choosing the more luminous stars first. Note that 
the final list is not the 150 brightest, because these are 
often spatially clustered and do not span the full area. 



Next, we visually identify a roughly linear shift be- 
tween the two coordinate systems and apply the trans- 
formation to our subset of matching stars. This com- 
parison acts both as a visual check that the star lists we 
are using are well constructed, and also as a first pass at 
determining the final transformation. 

The final transformation is determined iteratively 
with the routine MATCH developed by Michael Rich- 
mond, based on th e method of triangles described in 
Va ldes et al.l (I1995D . First we find a linear fit between 
the two coordinate systems. We then use that solution 
as a starting point for a quadratic solution. We find that 
a cubic solution is generally unnecessary for the trans- 
formation between the distortion corrected WFC3 and 
ACS images. 

After determining the transformation between the two 
coordinate systems with our set of 150 matching stars, 
we apply the transformation to the entire NIR dataset, 
bringing it into the optical coordinate system. The final 
step is to then use a separation criteria to determine if 
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Figure 2. continued 

there is a good match. We find that a separation of 
0.07" works well across the entire field. Typically 90% 
of the stars in the NIR catalogue are well-matched to a 
star in the optical catalogue. Of the remaining 10%, the 
bulk are either located in the wings of saturated stars 
or in the chip gap in the ACS camera. The F160W vs. 
{F8UW - F160W) CMDs for all of the program galaxies 
are shown in Figure [2j 

3. THE CONTRIBUTION OF TP-AGB AND RHEB STARS 
TO THE NIR FLUX OF GALAXIES 

The primary goal of this paper is to constrain the con- 
tribution of late stage stellar evolution to the NIR fluxes 
of galaxies with different metallicities and star formation 
histories. To do so, we must (1) determine the total flux 
falling in the WFC3 F160W fiher for each galaxy, and (2) 
determine the flux from the TP-AGB and RHeB stars in 
the same area. 

Determining the total fluxes of the sample fields is ac- 
tually non-trivial, as many of the images do not contain a 



clean measure of the background sky. Although we can 
easily determine the flux of the individual bright stars 
in each frame, we cannot directly photometer the stars 
that are too faint to detect. Instead we chose to model 
the contribution from the faint end, by generating syn- 
thetic NIR CMDs based on the optically derived SFHs, 
as described below. We will use these synthetic CMDs 
to both construct a total flux for each field, and to test 
model prescriptions for the TP-AGB and RHeB stars. 

3.1. Synthetic NIR CMDs 

We create synthetic NIR CMDs for each galaxy based 
on the SFHs derived from very deep optical HST imaging 
(see Section [2TT|) . We input the measured SFHs, redden- 
ing values, and distance moduli into CalcSFH to produce 
model Hess diagrams in the F81AW and F160M^ filters. 
We then sample th ese model He ss Diagrams with the 
routine NoisyCMD ()Dolphinll2002[ ) to generate synthetic 
photometry of each galaxy field down to K dwarfs^ 

NoisyCMD uses the Padova isochrones ((Marigo et al.l 
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120081 ) with updated bolometric corrections and Teff-color 
relations (jGirardi et alJ 120081 ) to populate the model 
CMDs. However, NoisyCMD does not include the ef- 
fects of long; lived thermal-pulses (i.e. 10,000 years, see 
iWagenhuber k, Groenewegeiilll998D which can scatter up 
to 20% of the lower mass TP-AGB stars to lower lumi- 
nosities. 

The Pad qva isochrones have been discussed in detail 
previously (iGirardi et all 120001 iMarigo fc Girardil 120071: 
iMarigo et all I2008D : here we include a brief descrip- 
tion. The primary distinctio n of the Padova isochrones 
over previous efforts (e.g.. iBruzual fc Charlod [200l 
iRaimondo et al.l 120051 : 1 Vazquez fc Leitheredl2005| ) is the 



detailed characterization of several key aspects of the TP- 
AGB phase, including: hot bottom burning, third dredge 
up, and variable atmospheric opacities. These effects are 
crucial for tracking the evolution of TP-AGB stars, es- 
pecially across the transition from oxygen rich to carbon 
rich phases, and the production of circumstellar dust. 
In addition, the Padova isochrones incorporate mass-loss 



from dust driven winds (|Winters et al.ll2000l [200l . and 
follow the TP-AGB evolution through the loss of the 
outer gaseous envelope. Circumstellar dust brings its 
own complications, and different pr escriptions can lead to 
different outcomes as explained in IMarigo et al.l (pOOS) . 
In the present work, we adopt the isochrones without 
circumstellar dust, however, briefly discussing the ef- 
fects dust may cause in the star counts and integrated 
fluxes. The bolo metric corrections are described in detail 
in Girardi et al . (2002, 2008). They were generat ed from 
the sp ectrophotometric standards assembled in iBohlinI 
(12007 ), and a large lib rary of spectral fluxes as sembled 
in IGirardi etHI (|2002| ) and lAringer et al.l (|2008( ). 

In an effort to account for any systematic offset be- 
tween the models and the data, we produce two itera- 
tions of the model photometry. In the first iteration, we 
generate model photometry of the more luminous stars. 
We compare the numbers of luminous RGB stars in the 
model to the number of RGB stars in the data (see brown 
box in Figurc[2]), and calculate a scaling between the two. 
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Figure 2. continued 

Wc then re-run the models to very faint levels, applying 
the scaling to the SFHs. This assures a good match be- 
tween the models and the data on the RGB. The scalings 
we calculate are typically less than 10% and give some 
indication of the uncertainty on the total luminosity we 
are measuring for each galaxy. 

We create two different synthetic CMDs. The first 
(middle pa nel of Figure El) is b ased on the 2008 Padova 
isochrones ()Marigo et al" [20081) , which form the basis for 
several commonly adopted stellar population synthesis 
codes in use today. The s econd includes the u pdated 
Padova TP-AGB models of IGirardi et"all ([Mot ). which 
has a shorter lifetime for low mass, low metallicity TP- 
AGB stars. This latter model effectively lowers the to- 
tal number of TP-AGB stars in the model CMDs (right 
panel of Figure [2]). 

The model CMDs are shown in Figure [2] and (by de- 
sign) are well-matched to the observed CMDs at the RGB 
(see Table 2). 



3.2. Total 1.6 iim Fluxes 

While the model CMDs do a good job reproducing 
the well understood main sequence and red giant branch 
stars, they do not necessarily reproduce the most lumi- 
nous stars which are in phases of late stage stellar evo- 
lution (i.e. the TP-AGB and He burning phases that we 
are investigating). We therefore adopt a total flux for 
each galaxy based on a hybrid of model -|- data fluxes, 
such that: 

Total Stellar Flux — Observed Flux of Luminous Stars -f 
Model Flux of Faint Stars (1) 

We chose a splice point between the data and the model 
at F160M^ — 23 mags, typically over a magnitude below 
the TRGB in F160W. At this flux level, the HST mag- 
nitudes are well measured, and are good to within 0.05 
mags. This splice point is also much brighter than the 
typical 50% completeness limit of F160W « 25 mags. 
When calculating the observed flux, we exclude the 
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Table 2 

Total Flux and Numbers of RGB Stars in Each Galaxy 



Galaxy 


Total F160W Flux 


# RGB 


# RGB 


# RGB 




ergs cm 


2 s-1 


data 


model 2008 


model 2010 


DD071 


7.550e- 


16 


1958 ± 44 


1940 ± 44 


1906 ± 43 


DD078 


1.330e- 


15 


3050 ± 55 


2913 ± 53 


3166 ± 56 


DD082 


4.570e- 


15 


9605 ± 98 


9320 ± 96 


9581 ± 97 


ESO540-030 


4.090e- 


16 


741 ± 27 


756 ± 27 


728 ± 26 


HS117 


3.980e- 


16 


844 ± 29 


805 ± 28 


899 ± 29 


IC2574-SGS 


4.570e- 


15 


7187 ± 84 


7299 ± 85 


7018 ± 83 


KDG73 


l.SlOe- 


16 


449 ± 21 


452 ± 21 


370 ± 19 


KKH37 


6.830e- 


16 


1214 ± 34 


1176 ± 34 


1293 ± 35 


M81-DEEP 


9.590e- 


16 


1032 ± 32 


1024 ± 32 


1076 ± 32 


NGC0300-WIDE1 


6.390e- 


15 


2516 ± 50 


2409 ± 49 


2516 ± 50 


NGC2403-HALO-6 


8.190e- 


16 


695 ± 26 


716 ± 26 


679 ± 26 


NGC2976-DEEP 


1.450e- 


15 


2879 ± 53 


2902 ± 53 


2916 ± 54 


NGC3077-PHOENIX 


1.270e- 


15 


2116 ± 46 


2156 ± 46 


2176 ± 46 


NGC3741 


l.lOOe- 


15 


1514 ± 38 


1524 ± 39 


1490 ± 38 


NGC4163 


3.780e- 


15 


5088 ± 71 


4981 ± 70 


5140 ± 71 


NGC7793-HALO-6 


6.920e- 


16 


1033 ± 32 


1059 ± 32 


1017 ± 31 


SCL-DEl 


3.000e- 


16 


588 ± 24 


548 ± 23 


629 ± 25 


UGC4305-1 


3.410e- 


15 


3664 ± 60 


3588 ± 59 


3725 ± 61 


UGC4305-2 


3.250e- 


15 


4087 ± 63 


4072 ± 63 


4064 ± 63 


UGC4459 


1.250e- 


15 


1815 ± 42 


1771 ± 42 


1787 ± 42 


UGC5139 


1.500e- 


15 


2783 ± 52 


2821 ± 53 


2647 ± 51 


UGC8508 


1.960e- 


15 


2094 ± 45 


2053 ± 45 


2097 ± 45 


UGCA292 


2.470e- 


16 


354 ± 18 


324 ± 18 


341 ± 18 



^ Data brighter than F160W= 23 mag plus model fainter than F160W= 23. 
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small number (if any) of extremely IR bright stars, 6 — 8 
magnitudes brighter than the TRGB (set by eye for this 
paper). These stars are likely to be foreground stars 
(see Section 3.6 below). Their flux is also excluded 
from the calculated contribution from the TP-AGB and 
RHcB phases. In general the foreground numbers are 
expected to be small in the TP-AGB and RHeB regions 
of the CMP for f ields of this small angular size (e.g. 
IGirardi et afllMol) . 

To give an example of the type of calculation we will 
be making in the following sections, Figure [3] shows the 
stellar (model + data) cumulative F160W flux fraction 
as a function of stellar magnitude for program galaxy 
UGC4305-1. This plot is divided into three regions. The 
left-most region shows the contribution of stars brighter 
than the TRGB. The right-most region gives the con- 
tribution from faint MS stars. The central region is 
dominated by RGB stars, but also contains some MS, 
fainter AGB, and fainter core helium burning stars. This 
plot demonstrates one of the reasons why galaxy model- 
ers have preferred to use NIR fluxes to constrain stellar 
masses. The bulk of the light is from well-modeled RGB 
and MS stars. However, even in this galaxy, which has a 
well-developed RGB, stars brighter than the TRGB con- 
tribute more than 30% of the light. At high-z— where the 
RGB has had little time to develop, we expect that the 
NIR luminous RHeB and TP-AGB stars will contribute 
significantly larger fractions to the total. 

The total F160Vt^ fluxes for each galaxy are given in 
Table El 

3.3. Identifying Late Stage Stellar Evolution Sequences 

We identify stars on the TP-AGB and RHeB sequences 
by selecting them in color-magnitude space. For the pur- 
poses of this paper, we are interested in those stars that 
might affect stellar mass estimates of high redshift galax- 
ies, and thus focus on stars more luminous than the 
TRGB. AGB stars fainter than these limits exist, but 
will not contribute significantly to the total light because 
their numbers will be dwarfed by the more numerous 
longer-lived RGB stars. 

Ideally we would use the IR observations to select a 
complete sample of the NIR luminous TP-AGB stars in 
each system. The TP-AGB stars are more luminous in 
the IR bands compared with the optical data, and their 
colors are less affected by dust reddening than in the 
optical. Thus a more complete census should be possible 
in the NIR compared to the optical, which may actuall y 
miss large numbers of TP-AGB stars (|Bover et al.ll2009D . 
Note: even the NIR can mis s the most dust obscured 
AGB stars fBover et al.l 120091) . but these stars will not 
affect the NIR luminosities of our galaxies because they 
are NIR faint. 

Unfortunately, there is a problem with using the NIR 
CMDs to cleanly identify TP-AGB stars; the TP-AGB 
sequence above the TRGB has a similar IR color as the 
younger RHeB sequence thus making these two stellar 
classes hard to distinguish (Figure[T]). To cleanly identify 
a complete set of luminous TP-AGB stars, we therefore 
select them in the {F8UW - F160W) CMDs (Figure[2]), 
which have a much larger color separation between the 
TP-AGB and RHeB sequences. 

To select TP-AGB and RHeB stars, we create regions 
in the CMD space that isolate these sequences (Figure[2]). 



To define the regions, we use a galaxy which has well pop- 
ulated TP-AGB and RHeB sequences, UGC 4305 (Figure 
[T]). As discussed above, the faint-end limit is set by the 
TRGB, and the bright limit is set to exclude luminous 
foreground stars. We shift these regions for each sub- 
sequent galaxy, applying a vertical shift to account for 
differences in distance and TRGB flux, and a horizon- 
tal shift to account for reddening variations. In addition 
to the RHeB and TP-AGB sequences, we also include a 
box that contains a large fraction of the upper RGB se- 
quence. Tables [2] -in give the numbers of luminous RGB, 
AGB, and RHeB stars in each galaxy. 

We find that the {F8UW - F160W^) CMDs contain 
roughly the same number of luminous stars (brighter 
than the TRGB) as the NIR only CMDs. Thus, we are 
unlikely to be missing large numbers of TP-AGB stars 
in the final analysis, although rare highly reddened TP- 
AGB stars could be absent. 

3.4. Fraction of NIR Light Produced by TP-AGB Stars 

Figure |4] shows the fraction of the F160W light pro- 
duced by the TP-AGB as a function of the fraction of 
young stars in each galaxy (as estimated from the SFH 
routine CalcSFH). The fractional flux contribution from 
TP-AGB stars in this sample varies from ~ 1% to 17%, 
with a trend of increasing contribution by the TP-AGB 
with an increasing fraction of young stars. Uncertain- 
ties on the flux fractions are derived from the Poisson 
uncertainties of the numbers of TP-AGB stars and their 
typical fluxes. 

Figure HI considers recent star formation on two 
timescales — 2 Gyr (left), and 0.3 Gyr (right). The 
shorter timescale tracks the lifetime of the most mas- 
sive TP-AGB stars (e.g. M > 3 Mq), while the longer 
timescale tracks the more common but less massive (e.g. 
M = 2 - 3 Mq) TP-AGB stars. Interestingly the scat- 
ter in the plot is reduced for the shorter star formation 
timescale. The reduced scatter in this second version of 
the plot may be indicating that the trend is driven by the 
most massive TP-AGB stars. Alternatively it may just 
be showing that the SFHs are better constrained for the 
youngest ages. Larger samples with better constrained 
SFHs at older ages could be useful for explaining this 
interesting result. 

While there is a strong trend with fractional age, there 
does not appear to be an equivalent trend with metal- 
licity. The metallicity shown in this plot is the mean 
expected metallicity for stars that are 1 Gyr old (as es- 
timated by CalcSFH). Low and high metallicity systems 
both appear to be following the same general trends of 
increasing TP-AGB contribution with increasing fraction 
of young stars. 

There is evidence that SFH uncertainties are contribut- 
ing to the scatter in the left hand version of this plot. 
Two galaxies, ESO540-030 and HS117, in particular ap- 
pear discrepant in the left-hand panel of Figure |31 show- 
ing less TP-AGB light than their SFH might imply. How- 
ever, the SFH is highly uncertain for ESO540-030. For 
instance the best fit SFH predicts a moderately high 
metallicity ([m/H] > —0.8) for the youngest populations 
of this low-mass galaxy. This is likely an over-estimate, 
as the metallicity is significantly lower for the bulk of cos- 
mic time. While other studies have also suggested higher 
metallicities for the most recent stars (Jerien &: Rcikubal 
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Figure 2. continued 

120011 ). the values they derive are still [m/H] < — 1. In 
the case of HS117 a handful of extremely bright RHeB 
stars may also be lowering the AGB contribution to the 
total. As we discuss in Section 3.6, these are likely fore- 
ground stars, further complicating the interpretation of 
these results. However, when we replot this Figure, now 
against the fraction of stars younger than 0.3 Gyrs (right 
hand panel of Figure |4|, ESO540-030 and HS117 are no 
longer significantly deviant. 

Another galaxy, UGCA292, is also somewhat anoma- 
lous in Figure ID UGCA929 is predicted to have the 
largest fraction of young stars, but its TP- AGB popula- 
tion does not account for a correspondingly large flux 
fraction compared with the other galaxies. However, 
UGCA292 is the least-populated galaxy in the sample 
(see Fig. 2), which causes not only the increased uncer- 
tainties depicted in the figure, but also potential difficulty 
in the derivation of the SFH. 

The observed TP- AGB fiux fractions are summarized 
in Table 131 Poisson uncertainties are also quoted. 



3.5. Fraction of NIR Light Produced by RHeB Stars 

Figure [5] shows the contribution of RHeB stars to the 
1.6 /im fiuxes of galaxies as a function of the fraction of 
young stars. Here we only plot the smaller age range, 
using the mass fraction of stars younger than 0.3 Gyrs. 
This represents the timescale (or progenitor mass range 
> 3.5 M0) over which RHcBs contribute significant frac- 
tions of the IR luminosity of galaxies. While less massive 
stars will go through a core helium burning phase, they 
will not reach luminosities brighter than the TRGB, and 
will instead populate the horizontal branch or red clump 
of the CMD. 

Again we see a similar trend where galaxies with a 
higher fraction of ongoing or recent star formation tend 
to show a larger contribution from the RHeB phases of 
stellar evolution. Interestingly, the contribution of RHeB 
stars can match or even exceed the contribution from the 
TP- AGB phase of stellar evolution, reaching as high as 
21% of the total in NGC 2403. While there is not a 
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Figure 2. continued 

strong trend in the RHeB flux fraction with metaUicity, 
there may be some favoring of lower flux fractions for 
galaxies with low metaUicity. Such a trend could indi- 
cate that RHeBs are rarer in low metaUicity systems. 
Comparing the red to blue helium burning fractions as 
a function of galaxy metaUicity could shed more light on 
this issue (|McQuinn et al.l[2011t ). However, significantly 
more galaxies, especially galaxies with larger fractions of 
young stellar populations, should be used to determine 
if this trend indeed exists. 

As with the TP- AGE stars HS117 is somewhat deviant 
in this plot (although now in the opposite direction as in 
the TP-AGB plot. In this handful (4 stars) of ex- 

tremely bright stars are pushing the RHeB flux fraction 
higher than 10% even though the SFH suggests little star 
formation at these young ages. With these small numbers 
it is hard to draw significant conclusions, as foreground 
stars could be important. In the following section we will 
attempt to statistically account for any foreground stars. 

The other galaxy that might be considered somewhat 



deviant in this plot is UGC 292. However, once again, 
the large uncertainty in the SFH at young ages mean 
that this galaxy could actually belong closer to the main 
locus of points. 
These results are also summarized in Table |4l 

3.6. Foreground Stars 

With any CMD studies, foreground (or background) 
contamination can pose a problem. There may be Milky 
Way stars that have similar colors and luminosities as 
the stars in the program galaxies. This effect could tend 
to artificially increase the numbers of stars on the RHeBs 
or TP-AGB sequences. We actually benefit from the fact 
that our fields are spatially very small, only 4.7 square 
arcmin. This area is much smaller than is needed to 
image a typical Local Group dwarf galaxy (> 200 square 
arcmin) and thus has much lower field contamination. 
However, contamination may still be important for the 
RHeB region, where the observed numbers can be small. 

To statistically estimate the foreground contamination. 
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Table 3 

AGB Star Properties of Each Galaxy 



Galaxy I # AGB fAGB/ftot 

data model 2008 model 2010 data model 2008 model 2010 



DD071 


146 ± 12 


3 


62 




1 


94 


1 


53 







82 





08 







01 


4 


48 


± 


2 


99 


1 


96 


± 


1 


31 


DD078 


273 ± 16 


3 


14 


± 


1 


51 


2 


64 




1 


27 





08 







01 


2 


99 


± 


1 


68 


2 


81 


± 


1 


58 


DD082 


1046 ± 32 


2 


44 


± 


1 


33 


1 


07 


± 





59 





10 


± 





00 


2 


96 


± 


1 


35 


2 


22 


± 


1 


02 


ESO540-030 


69 


± 


8 


5 


86 


± 


2 


12 


4 


14 




1 


50 





07 


± 





01 


8 


51 


± 


4 


33 


6 


70 


± 


3 


41 


HS117 


63 ± 7 


2 


57 


± 


1 


09 


2 


05 


± 





87 





06 


± 





01 


3 


20 


± 


1 


44 


2 


85 


± 


1 


28 


IC2574-SGS 


1504 ± 38 


1 


72 







63 


1 


03 







38 





17 


± 





00 


2 


22 


± 





59 


1 


97 


± 





52 


KDG73 


60 


± 


7 


1 


82 







69 


1 


25 







47 





13 


± 





02 


2 


13 


± 


1 


21 


1 


49 


± 





85 


KKH37 


109 


± 


10 


2 


62 




1 


14 


1 


34 







58 





10 


± 





01 


2 


32 


± 


1 


04 


1 


54 


± 





69 


M81-DEEP 


157 


± 


12 


4 


79 







73 


4 


94 







75 





08 


± 





01 


5 


89 


± 





93 


7 


43 


± 


1 


17 


NGC0300-WIDE1 


465 


± 


21 


1 


99 







51 


1 


80 







46 





15 


± 





00 


2 


61 


± 





58 


2 


75 


± 





61 


NGC2403-HALO-6 


182 


± 


13 


1 


69 


± 





71 


1 


37 


± 





58 





17 


± 





01 


1 


74 


± 





70 


1 


80 


± 





72 


NGC2976-DEEP 


293 


± 


17 


2 


76 


± 


2 


01 


1 


27 







92 





10 


± 





01 


3 


00 


± 


2 


49 


1 


43 


± 


1 


19 


NGC3077-PHOENIX 


215 


± 


14 


2 


91 


± 


1 


20 


2 


65 




1 


09 





09 


± 





01 


2 


74 


± 


1 


15 


2 


65 


± 


1 


11 


NGC3741 


233 


± 


15 


2 


07 


± 





54 


1 


30 


± 





34 





15 


± 





01 


2 


60 


± 





82 


2 


04 


± 





64 


NGC4163 


640 


± 


25 


2 


23 




1 


26 


1 


06 







60 





12 


± 





00 


2 


48 


± 


1 


42 


1 


67 


± 





96 


NGC7793-HALO-6 


150 


± 


12 


2 


09 


± 





87 


1 


12 


± 





47 





15 


± 





01 


2 


38 


± 





90 


1 


75 


± 





66 


SCL-DEl 


67 


± 


8 


1 


67 


± 


1 


14 





94 


± 





64 





07 


± 





01 


1 


53 


± 


1 


19 


1 


11 


± 





87 


UGC4305-1 


740 


± 


27 


1 


94 







44 


1 


50 







34 





15 


± 





00 


2 


81 


± 





66 


3 


17 


± 





75 


UGC4305-2 


721 


± 


26 


2 


11 


± 





58 


1 


47 







40 





17 


± 





00 


2 


75 


± 





69 


2 


64 


± 





66 


UGC4459 


262 


± 


16 


3 


14 







77 


2 


53 







62 





13 


± 





01 


2 


82 


± 





90 


2 


82 


± 





90 


UGC5139 


439 


± 


20 


2 


43 




1 


05 


1 


39 







60 





14 


± 





01 


3 


21 


± 


1 


20 


2 


38 


± 





89 


UGC8508 


268 


± 


16 


2 


37 


± 





91 


1 


24 


± 





48 





14 


± 





00 


2 


70 


± 


1 


21 


1 


67 


± 





75 


UGCA292 


62 


± 


7 


2 


81 


± 





77 


2 


00 


± 





55 





17 


± 





02 


3 


69 


± 


1 


06 


4 


39 


± 


1 


26 



Table 4 

RHeB Star Properties of Each Galaxy 



Galaxy I # RHeB fRH.B/ftot ^-ff 

data model 2008 model 2010 data model 2008 model 2010 



DD071 


2 ± 1 





50 


± 


1 


87 


0.50 




1.87 





00 


± 


0.00 





00 


± 


0.00 





43 


± 


0.86 


DD078 


5 ± 2 





20 


± 





81 


0.20 




0.81 





00 


± 


0.00 





00 


± 


0.00 





44 


± 


0.88 


DD082 


39 ± 6 


1 


51 


± 


2 


74 


0.49 




0.88 





02 


± 


0.00 





22 


± 


0.33 





23 


± 


0.19 


ESO540-030 


10 ± 3 





10 


± 





09 


0.10 




0.09 





03 


± 


0.01 





00 


± 


0.00 





06 


± 


0.12 


HS117 


8 ± 2 





12 


± 





12 


0.25 




0.23 





11 


± 


0.01 





01 


± 


0.00 





02 


± 


0.03 


IC2574-SGS 


479 ± 21 





46 


± 





13 


0.43 


± 


0.13 





19 


± 


0.01 





30 


± 


0.08 





23 


± 


0.08 


KDG73 


18 ± 4 





17 


± 





07 


0.44 


± 


0.19 





07 


± 


0.04 





13 


± 


0.14 





36 


± 


0.29 


KKH37 


15 ± 3 





40 


± 





24 


0.27 


± 


0.16 





05 


± 


0.01 





14 


± 


0.21 





12 


± 


0.10 


M81-DEEP 


16 ± 4 





38 


± 





90 


0.06 


± 


0.15 





02 


± 


0.01 





33 


± 


0.66 





30 


± 


0.30 


NGC0300-WIDE1 


117 ± 10 





44 


± 





13 


0.38 




0.12 





09 


± 


0.00 





33 


± 


0.24 





43 


± 


0.25 


NGC2403-HALO-6 


50 ± 7 





34 


± 





17 


0.42 


± 


0.21 





21 


± 


0.02 





24 


± 


0.15 





20 


± 


0.12 


NGC2976-DEEP 


1 ± 1 


8 


00 


± 


4 


79 


20.00 


± 


11.96 





00 


± 


0.00 


21 


88 


± 


11.54 


19 


80 


± 


14.85 


NGC3077-PHOENIX 


44 ± 6 





50 


± 





17 


0.59 




0.20 





05 


± 


0.01 





65 


± 


0.26 





72 


± 


0.29 


NGC3741 


63 ± 7 





71 


± 





18 


0.68 


± 


0.18 





09 


± 


0.01 





46 


± 


0.14 





41 


± 


0.15 


NGC4163 


120 ± 10 





27 


± 





12 


0.41 


± 


0.18 





08 


± 


0.01 





22 


± 


0.12 





23 


± 


0.09 


NGC7793-HALO-6 


29 ± 5 





52 


± 





17 


0.62 


± 


0.21 





10 


± 


0.01 





14 


± 


0.19 





10 


± 


0.20 


SCL-DEl 


5 ± 2 





20 


± 





19 


0.20 


± 


0.19 





01 


± 


0.01 





00 


± 


0.00 





08 


± 


0.08 


UGC4305-1 


296 ± 17 





68 


± 





15 


0.64 


± 


0.14 





16 


± 


0.01 





45 


± 


0.10 





41 


± 


0.12 


UGC4305-2 


280 ± 16 





67 


± 





14 


0.73 


± 


0.16 





10 


± 


0.01 





66 


± 


0.15 





56 


± 


0.17 


UGC4459 


49 ± 7 


1 


08 


± 





33 


1.08 


± 


0.33 





10 


± 


0.01 





42 


± 


0.18 





48 


± 


0.21 


UGC5139 


167 ± 12 





68 


± 





21 


0.68 




0.21 





09 


± 


0.02 





71 


± 


0.32 





82 


± 


0.28 


UGC8508 


49 ± 7 





63 


± 





37 


0.41 




0.24 





05 


± 


0.01 





59 


± 


0.40 





75 


± 


0.26 


UGCA292 


22 ± 4 


1 


18 


± 





40 


1.23 


± 


0.41 





16 


± 


0.03 





47 


± 


0.17 





59 


± 


0.22 



we run TRILEGAL (jGirardi et alJ l2005f ) which mod- 
els the Milky Way contamination for a given field size, 
in a given input direction. We run TRILEGAL with 
the canonical settings including a thin disk component, 
a bulge component, and a halo component. We only 
consider model foreground stars that are brighter than 
F160M^ = 23 mag, roughly the magnitude cutoff for real 
data in our calculation of the total fluxes of the program 
galaxies. For every star in the foreground model, we 
determine if there is a real star within 0.5 mags in CMD- 
space. If there is, we flag the closest one as a potential 



foreground star. Thus we only account for plausible fore- 
ground stars. For instance a model foreground star that 
is 0.5 mags brighter than the brightest actual star will 
be assigned the flux of the actual star. Likewise, model 
stars that are brighter still, will not be considered at all. 

Table [5] presents statistical estimates for the numbers 
and fluxes of foreground stars in the direction of each 
galaxy in our sample. The foreground stars have been 
classified by the region of the CMD in which they are 
found, i.e. RHeB and TP- AGB. An estimate of the total 
flux from all plausible foreground stars is also provided 
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UGC4305-1 




F1 60W [mag] 

Figure 3. The cumulative F160W flux fraction as a function of 
stellar magnitude for program galaxy UGC4305-1. The contribu- 
tion from stars brighter than the TRGB is given to the left of the 
first vertical line and includes light from TP-AGB and RHeB stars. 
The contribution from the faint MS is given to the right of the sec- 
ond vertical line. The central region is dominated by RGB stars 
but also contains contributions from brighter MS, fainter AGB and 
fainter core helium burning stars. The bulk of the light is from 
well modeled MS and RGB stars, however, even in this galaxy 
more than 30% of the F160W light is from stars brighter than the 
TRGB, namely TP-AGB and RHeB stars. 

in the table, including stars that are not in the RHeB 
and TP-AGB regions. 

The TRILEGAL models typically predict fewer than 5 
plausible foreground stars in the TP-AGB region of the 
CMD. Although for KKH 37 there may be twice that 
number. Foreground stars account for fewer than 5% of 
the AGB flux in these galaxies, except for KKH 37 where 
the foreground may be contributing as much as 20%. 
Because the foreground contribution is smaller than the 
Poisson uncertainties on the numbers of TP-AGB stars, 
we will not make any special attempt to account for it in 
our additional analysis of the TP-AGB. 

Typically there are on the order of 5 foreground stars 
in the RHeB region, with as many as 10 in NGC 2403. 
As half of the program galaxies contain fewer than 30 
stars in the RHeB region, foreground could potentially 
be responsible for a large fraction of the measured RHeB 
flux. For instance, HS117, which is thought to have little 
recent star formation, is found to have 8 RHeB stars 4 
of which are very luminous and make up the over 90% of 
the flux. However, the TRILEGAL model predicts that 
at least three of those luminous stars are foreground. 
The one remaining RHeB star could also well be fore- 
ground given the small number statistics. Removing the 
foreground stars from the RHeB sequences actually im- 
proves the correlation between the flux contribution of 
RHeBs and the star formation histories of the galaxies 
as shown in the right hand panel of Figure [5j After fore- 
ground correction HS117 is no longer a deviant point 
in the plot. Because of the uncertainties introduced by 
foreground stars on the RHeB sequences, we will only 
consider galaxies with more than 30 RHeBs for the re- 
mainder of the RHeB analysis. 

3.7. Comparison with Models 



In addition to quantifying the flux contributions of TP- 
AGB and RHeB stars to the NIR luminosities of galaxies, 
we check if these results can be reproduced by SPS codes 
based on the 2008 and 2010 Padova isochrones. For each 
galaxy, we use the best-fit optically derived SFH to model 
the stellar content in the NIR. We then compare the 
numbers and fluxes of the synthetic TP-AGB and RHeBs 
to the real data. 

To determine the uncertainties on the numbers and 
fluxes of synthetic TP-AGB and RHeB stars, we create 
100 additional model CMDs for each galaxy. Each of 
these 100 models is created with a different SFH chosen 
to span the range of acceptable SFHs as determined by 
CalcSFH. The adopted uncertainty of the model mea- 
surement is then given by the standard deviation of the 
100 Monte Carlo results. 

Figure [6] shows an example model/data comparison for 
UGC 4305-1. The four panels plot the observed numbers 
(top) and flux contributions (bottom) of TP-AGB and 
RHeB stars for actual data (filled circles), and the best- 
fit SPS models from the 2008 (open circles) and 2010 
(diamonds) Padova isochrones. Also shown are the re- 
sults for 100 MC simulations based on the 2010 models 
(histograms) that span the full range of acceptable SFHs. 
For UGC4305-1 the SPS models tend to over-predict the 
observed numbers and fluxes of the TP-AGB stars and 
under-predict the observed numbers and flux contribu- 
tions of RHeB stars, although the uncertainties on the 
model results (i.e. the widths of the histograms) are 
large. 

Figure [7] compares the predicted-to-observed numbers 
of TP-AGB stars for the full sample of galaxies. Results 
from the best-fit SPS models based on both the 2008 (cir- 
cles) and 2010 (diamonds) Padova isochrones are shown. 
It is clear that 2008 models tend to overestimate the 
numbers of TP-AGB stars, presenting excesses in the 
numbers of TP-AGB stars by factors of 1.7 — 6. The 
situation is largely remedied in the 2010 models: indeed, 
for 65 % of the galaxies the ratio between modeled and 
observed numbers is consistent with unity to within the 
measured uncertainties. However, there remains a small 
systematic bias to larger numbers in the 2010 models by 
a weighted mean factor of 1.5 with a standard deviation 
0.5. The weighted mean factor drops to 1.4 ±0.5 if we re- 
move the two most discrepant galaxies. Considering the 
relatively large uncertainties in the SFH and metallicities 
of these galaxies, this agreement can be considered quite 
good. 

Figure [7] also compares the observed vs. predicted flux 
contribution of TP-AGB stars. The best-fit SPS models 
show a larger offset in TP-AGB flux than in numbers, 
with models generally over-predicting the flux contribu- 
tion. The typical offset for the 2008 models is a factor 
three, but as expected, the 2010 models do a better job 
of reproducing the data, especially for galaxies with lit- 
tle on-going star formation. However, for galaxies with 
recent star formation, the 2010 models appear to be just 
as discrepant as the 2008 models, with offsets reaching a 
factor of 2 - 3 or larger. The weighted mean offset be- 
tween the model/data flux ratio is 2.3 with a standard 
deviation of 0.8. 

Two galaxies present very high model/data discrepan- 
cies for both number and flux contribution of AGB stars, 
ESO540-030 with a model/data flux ratio = 6.70 and 
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Figure 4. The fraction of F160W flux produced by TP-AGB stars as a function of the young stellar populations in each galaxy. Mass 
fractions of stars less than 2 Gyrs (typical age range for TP-AGB brighter than TRGB) is shown on the left, while mass fractions of stars 
less than 0.3 Gyrs (age range for the most massive and luminous TP-AGB stars) is shown on the right. Galaxies with a higher fraction 
of young stars tend to show a larger contribution from TP-AGB stars then more evolved galaxies. The TP-AGB contribution reaches as 
high as 17% even in this sample which contain a well developed red giant branch. In high-2 galaxies where the RGB has not developed, 
the TP-AGB contribution is likely to be significantly higher. We examine the influence of metallicity on the trend, using the expected 
metallicity of the 1 Gyr old population from the measured SFHs, and find no obvious trend with metallicity. Galaxies ESO540-030 and 
HS117 are somewhat deviant in the first version of this plot but not in the second, demonstrating the difflculties in interpreting the effects 
of rare but luminous stars on a stellar population. 



M81-DEEP with model/data flux ratio = 7.43. These 
galaxies are discussed further in Section 4.1. The third 
most discrepant galaxy, UGCA292 with model/data flux 
ratio = 4.39, is the least populated galaxy in our sam- 
ple. Therefore its large measurement uncertainties are 
driven by both small number statistics and large uncer- 
tainties in the SFH. If we remove these galaxies from the 
weighted mean we derive an overall model/data flux ratio 
of 2.2 ±0.8. 

Figure [8] compares the data to the SPS model pre- 
dictions for RHeB stars. In this case, the best-fit SPS 
models tend to under-predict the numbers and flux con- 
tribution of stars; both are under-estimated by a mean 
factor of 2.0 ±0.6 (for galaxies with larger than 30 RHeB 
stars, shown as filled diamonds in the plot). There is no 
appreciable difference between the SPS models based on 
the 2008 and 2010 Padova isochrones, so only the 2010 
models are shown in this figure. Likewise, correcting for 
foreground contamination does not alter these results as 
they are based on the galaxies with the largest number of 
RHeB stars. Explanations for the model/data differences 
will be explored in the discussion section. 

4. DISCUSSION 

Our results show that both the TP-AGB and RHeB se- 
quences can contribute significantly to the NIR fiux of a 
galaxy. Even for galaxies with well-developed red giant 
branch populations, the combination of these two late 
stages of stellar evolution can make up almost 40% of 



the NIR light, while comprising negligible stellar mass. 
As a result, these phases must be well-calibrated to accu- 
rately estimate stellar masses of galaxies. Making mat- 
ters worse, the parent populations for RHeBs and the 
most luminous TP-AGB stars are comprised of massive 
(e.g. M > 3.5 M0) young stars with short fifetimes. We 
therefore expect significant variations in the NIR M/L 
ratio of galaxies on very short timescales (e.g. < 300 
Myrs), especially in the early universe, where these stars 
will dominate the light. 

In addition to affecting M/L ratios, these rare but lu- 
minous populations may also be responsible for some 
of the scatter in key NIR galaxy scaling relatio ns such 
as NIR metallicity-luniinosity relations (e.g. Salze r et al.l 
2005) , and NIR vers ions of the TuUy-Fisher relation (e.g. 
iConselice et al.|[2005l ) . Typically the scatter in these NIR 
relations is smaller by ~ 30% compared to the optical, 
from avoiding the large scatter induced by dust obscu- 
ration and very young blue stars. However, significant 
scatter remains — e. g. more than the f ormal uncer- 
tainty in metallicity (jSalzer et al.l 120051 ). iSalzer et al.l 
(|2005f ) demonstrate that the scatter in the NIR does not 
correlate with the instantaneous star formation rate (as 
measured by Balmer line strengths). However, some of 
the remaining scatter could well be attributed to differ- 
ences in the ratios of evolved luminous stars. For in- 
stance we show here that galaxies with recent star for- 
mation and large RHeB populations could be as much 
as 20% brighter at 1.6 fim compared with similar mass 
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Figure 5. Same as Figure |4] only now for tiie RHeB stars. Here we only plot against the shorter timescale of 0.3 Gyrs, as this is the 
timescale for RHeBs. The left panel shows the raw results, while the right applies a statistical correction for foreground contamination as 
estimated by TRILEGAL (see Table [Sjl. Again there is a correlation, RHeBs contribute a larger fraction of the FldOW flux in galaxies 
with young stellar populations. The flux contributions from RHeBs can be as large as those from TP-AGB stars and reach as high as 21% 
of the total, or 18% after a foreground correction. As with the TP-AGB there is no strong trend in RHeB flux fractions with metallicity 
although there may be a slight favoring of larger RHeB flux fractions for more metal rich systems. 



Table 5 

Statistical Estimates of Foreground Star Contamination from TRILEGAL 



Galaxy TP-AGB Region RHeB Region All Foreground^ 





# 


flux'' 


# 


flux'' 


# 


flux'' 


DD071 


3 


1.56e-18 


1 


1.36e-18 


9 


5.01e- 


18 


DD078 


3 


1.96e-18 


1 


9.44e-19 


6 


3.62e- 


18 


DD082 


2 


1.47e-18 


8 


4.81e-17 


14 


5.04e- 


17 


ESO540-030 


1 


6.91e-19 


3 


5.37e-18 


7 


6.52e- 


18 


HS117 


1 


1.98e-19 


4 


3.21e-17 


7 


4.66e- 


17 


IC2574-SGS 


4 


6.59e-18 


6 


5.35e-17 


15 


6.78e- 


17 


KDG73 


3 


1.12e-18 


6 


8.33e-18 


12 


1.37e- 


17 


KKH37 


11 


1.21e-17 


5 


2.64e-17 


25 


3.99e- 


17 


M81-DEEP 


3 


1.23e-18 


3 


4.17e-18 


9 


5.87e- 


18 


NGC0300-WIDE1 


2 


4.85e-18 


3 


2.64e-17 


10 


3.80e- 


17 


NGC2403-HALO-6 





0.00e-(-00 


10 


4.83e-17 


10 


4.83e- 


17 
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1 


3.39e-18 





O.OOe+00 


7 


7.33e- 


18 
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5 


1.40e-17 


7 


1.54c- 


17 


NGC3741 


1 


2.54e-18 


2 


1.79e-17 


5 


2.13e- 


17 


NGC4163 


1 


2.66e-18 


5 


4.83e-17 


12 


5.21e- 


17 


NGC7793-HALO-6 


1 


1.64e-18 


8 


1.53e-17 


13 


1.79e- 


17 


SCL-DEl 


4 


1.15e-18 





O.OOe+00 


8 
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18 
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6 
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5 


4.17e-17 


18 


5.11e- 


17 
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6 


5.46e-18 


5 


4.53e-17 


19 


6.39e- 


17 


UGC4459 


5 


5.91e-18 


5 


3.63e-17 


13 


4.52e- 


17 


UGC5139 


4 


2.26e-18 


5 


1.81e-17 


13 


2.10e- 


17 
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4 


3.17e-18 


2 


3.08e-18 


9 
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18 
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6 
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Figure 6. A comparison of the observed numbers and fluxes 
of TP-AGB and RHcB stars to model predictions for galaxy 
UGC4305-1. The observed data are shown as filled circles. The 
SPS model results based on the 2008 (open circles) and 2010 
(diamonds) Padova isochrones are also shown. In addition, we 
plot the results from 100 MC simulations of the 2010 model (his- 
tograms) that span the range of acceptable SFH as determined by 
CALCSFH. While the observations are generally offset from the 
simulation results, the scatter in the simulations is typically large. 
We take the standard deviations of the 100 MC simulations as a 
measure of the uncertainty in the model results. 

galaxies with httle recent star formation. Likewise the 
younger TP-AGB stars also can significantly impact the 
NIR fl ux of a galaxy on lo nger timescales. As the scatter 
in the iSalzer et al.l ()2005[ ) metallicity- luminosity relation 
is ^ 0.2 dex, luminous evolved stars could well be playing 
a significant role. 

Unfortunately, we have shown that the best-fit SPS 
models of our sample galaxies — based on the optically 
derived SFHs — currently have difficulties recovering the 
NIR flux contributions from RHeB and TP-AGB stars. 
The best-fit SPS models tend to over-predict the NIR 
fiuxes from TP-AGB stars and under-predict the NIR 
fiuxes from RHeB stars. This latter discrepancy is par- 
ticularly worrisome as the short-lived luminous RHeB 
stars are susceptible to rapid variations in the SFR. 
For both the RHeB and TP-AGB, several factors may 
be contributing to the model/data discrepancies. Some 
may be driven by difficulties in measuring accurate SFHs 
for these galaxies, especially on the short timescales for 
which these stars live. Difficulties with SPS modeling 
(converting SFH into a CMD), including oversimplifica- 
tions in the description of the Te//-luminosity varia- 
tions during thermal pulse cycles and in the description 
of extinction by circumstellar dust, may also play a role. 
Additionally, the stellar evolution codes on which the 
SPS models are based may require further updates for 
these difficult to model stages of evolution. We now ex- 
amine the model/data discrepancies in more detail and 
explore the implications for studies of galaxies in the lo- 
cal and high redshift universe. 

4.1. TP-AGB Model and Data Discrepancies 

iGirardi et all ()2010[ ) showed that SPS models, based 
on the 2008 Padova isochrones, over-predicted the num- 



bers of TP-AGB stars in optical CMDs of ten old, metal 
poor low-mass galaxies. They found that the predicted 
numbers of TP-AGB stars could be reconciled with the 
data if the lifetimes of the low-mass TP-AGB stars were 
significantly reduced from roughly 4 Gyrs to roughly 1 
Gyr, i.e. reduced to lifetimes similar to those of more 
massive AGBs. This conclusion formed the basis of the 
newly released Padova 2010 isochrones. 

Now we have expanded this investigation to 23 galax- 
ies with a wide range of SFHs, including dwarfs and spi- 
rals, and have examined th e behavior of t he TP -AGB 
in the NIR. We confirm the IGirardi et all (|2010( ) result 
that SPS models based on the 2008 Padova isochrones 
over- predict the numbers of TP-AGB stars. We also find 
that SPS models based on the 2010 version of the Padova 
isochrones have largely eliminated the over-prediction of 
TP-AGB numbers. While the new SPS models still show 
a small systematic bias to larger numbers than the data, 

65% of the sample galaxies have model number pre- 
dictions that overlap the data to within the uncertainties 
(Figure Ej). 

Figure |9| shows this comparison in more detail. For 
each galaxy, we plot histograms of the numbers of TP- 
AGB stars in the data and compare to the results from 
the 2008 and 2010 models, only now the model his- 
tograms are sub-divided by the age (mass) of the syn- 
thetic star. The 2010 model effectively eliminates the 
oldest (lowest mass) TP-AGB stars bringing the pre- 
dicted numbers of TP-AGB stars more in line with the 
data. We stress that it may be possible to reduce the 
2010 model numbers further while using the same 2010 
Padova isochrones, by means of reasonable changes to the 
SPS code: e.g., by adding temperature-luminosity varia- 
tions driven by the thermal pulse cycles, and/or dust ob- 
scuration variations. These effects alone could be able to 
reduce the model numbers by a good ~ 20% percent, and 
will be further explored by Rosenfield et al. (in prep.). 

We note that the most discrepant points, in the case of 
the 2010 models, correspond to the galaxies ESO540-030 
with an excess factor of 4.14, and M81-DEEP with 4.94. 
These two galaxies are discrepant in most of the plots 
presented in this paper. M81-DEEP is unusual in that 
it is the most metal-rich system in the sample. While 
the mean metallicity at 1 Gyr is listed at —0.4 in Ta- 
ble 1, this estimate is misleading. It hides the fact that 
most (75%) of the TP-AGB stars (which span a range of 
ages) have metallicities near solar (i.e. [AI/H] > —0.1), 
producing the very red TP-AGB branch in the modeled 
CMDs. At the very high metallicity end, it is possible 
that the models tuned to the LMC could have problems. 

Issues at high metallicity may also be affecting the 
modeled CMDs of several other galaxies in this sam- 
ple. For instance, galaxies DD082, IC2574-SGS, and 
NGC 300 all show plumes of very red TP-AGB stars 
in their modeled CMDs. These plumes are not nearly as 
well-populated or obvious in the observed CMDs of these 
galaxies. While the mean metallicities for the modeled 
TP-AGB populations in these galaxies are metal poor, 
these red plumes are metal rich (roughly solar). This 
result again suggests issues at the high metallicity end, 
either with the measured SFHs or the stellar evolution 
codes. 

The metalhcity of ESO540-030 is not unusually high, 
but its modeled CMD is characterized by a large fraction 
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Figure 7. A model/data comparison of the numbers and F160W fluxes of the TP-AGB stars in our sample galaxies. SPS model results 
based on the optically derived best-fit SFH and the 2008 Padova isochrones are shown as circles (left). SPS models based on the 2010 
Padova isochrones are shown as diamonds (right). Uncertainties on the model/data ratios based on 100 MC simulations that span the range 
of acceptable input SFHs are shown. While the 2008 models tend to over-predict the numbers of TP-AGB stars, especially for galaxies 
with little on-going star formation, the 2010 version of the models reproduce the data reasonably well; 65% of the models are equivalent to 
the data to within the measured uncertainties (see dotted 1/1 line). The weighted mean model/data number ratio (dashed line) is 1.5 (for 
the 2010 model) with a standard deviation of 0.5. While the updated best-fit SPS models do a reasonable job with the numbers of stars, 
they tend to over-predict the fluxes of the TP-AGB stars by a larger weighted mean factor of 2.3 with a standard deviation of 0.8, similar 
to the results from the 2008 model. 



of extremely red objects all along the RGB sequence, so 
the accuracy of its SFH may be called in question. If we 
neglect ESO540-030 and M81-DEEP, the overproduction 
problem of the 2008 models has been largely corrected in 
the 2010 models - even though the correction only applies 
to the old metal poor AGB stars. 

While we have largely accounted for the differences in 
TP-AGB number, there continues to be discrepancies be- 
tween the predicted and observed TP-AGB fluxes. How 
can we understand that the best-fit SPS models reason- 
ably reproduce the TP-AGB numbers, with typically less 
than ~ 50% excess, while the excess in F160W flux is 
much larger (> 200%)? The answer is that the excess 
occurs mainly in the form of relatively few but very lu- 
minous TP-AGB stars that are predicted to exist in the 
model but that are not observed in the data. This phe- 
nomenon can be seen in the CMDs of Figure 2. 



To better understand these predictions. Figure [TO] re- 
plots the 2010 models for 3 of the sample galaxies, now 
broken up into different time bins. The most luminous 
TP-AGB stars are found at the youngest ages, in the 
< 0.3 Gyr age bins, which correspond to the highest pro- 
genitor masses (M > 3.5 Mq). Additional luminous TP- 
AGB stars are found in the 0.3 — 1 Gyr age bins. Very few 
of the oldest (> 1 Gyr), or lowest mass, TP-AGB stars 
reach the brightness levels of their more massive counter- 
parts. The models appear to over-predict the numbers of 
the brightest TP-AGB stars (above the dashed line) by 
large factors, primarily at the youngest ages (< 1 Gyr). 

Unfortunately, due to a paucity of stars, constraints 
derived from TP-AGB stars in star clusters of the Mag- 
ellanic Clouds are relatively poor for the age i nterval 
between 0.1 and ~ 0.5 Gyr (jGirardi fc Marigol [20071. 
Therefore, Figures 2 and 1101 could be providing precious 
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Figure 8. Same as Figure[7] only now for the RHeB model/data comparison. Because there is no difference in the RHeB populations of 
the 2008 and 20f0 models only the 2010 version is shown here (diamonds). The best-fit SPS models tend to under-predict the numbers 
and fluxes of the RIfeB stars. The weighted mean fractional model/data differences for galaxies with significant numbers of RHeBs (> 30 
stars, flUed diamonds) is 0.5 ± 0.2, for both the numbers and fluxes of RHeBs. 



information at wtiere specifically the TP-AGB models 
need further improvement, and where additional obser- 
vational constraints are the most urgent. 

Other issues may also contribute significantly to the 
model/data mismatch. (1) Lack of a description for the 
pulse cycle luminosity variations (e.g. low-luminosity 
dip and flash luminosity) in the NoisyCMD code leads 
to a fraction (~ 20 — 30%) of the simulated AGB stars 
brighter than they should be. While the TP-AGB tracks 
in Marigo & Girardi (2007) do follow in detail the flash- 
driven luminosity variations, these features are not in- 
cluded in the stellar isochrones or in NoisyCMD, where 
the whole TP-AGB evolution is assigned the pre-flash 
maximum quiescent luminosity predicted by the core 
mass - luminosity relation. (2) Dust obscuration from 
circumstellar TP-AGB envelopes could affect both the 
observational selection criteria and/or the model predic- 
tions. In fact, the observational data could be "miss- 
ing" TP-AGB stars, with self-extinction hiding NIR lu- 
minous TP-AGB stars in the optical data (F814W) data 
we use to select them. Additionally, model prescriptions 



for the dustiest phases of the TP-AGB are not included 
in NoisyCMD, thus they may be predicted to be more 
luminous at NIR wavelengths, than if proper dust mod- 
eling were used. (3) The SFHs we have derived for our 
galaxies may be wrong and thus predict incorrect number 
of AGB stars. Accurate constraints on the SFHs on the 
short timescales of the most massive stars are difficult to 
obtain even from very deep CMDs. 

We now explore the last two possible limitations, and 
leave further examination of the stellar evolution codes 
to future papers in this series. 

4.1.1. Are the Fd,UW ~ Fl&QW CMDs Missmg NIR 
Luminous AGB Stars? 

iBoyer et al.l ((2009) demonstrated that optical searches 
will miss a large fraction of the most dust obscured AGB 
stars. As part of our search criteria, we have used some 
optical data {HST FSIAW), and therefore are likely 
to have missed the most dust obscured sources. How- 
ever, for our purposes we are only concerned about those 
missing TP-AGB stars that are actually luminous at 
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Figure 9. The numbers of TP-AGB stars in eacli galaxy (purple) 
compared to the 2008 (middle of each group of three) and 2010 
(right of each group of three) Padova model predictions. The total 
number of TP-AGB stars for each model are divided into differ- 
ent age bins, as indicated in the legend. In all cases the models 
over-predict the numbers of TP-AGB stars. The main difference 
between the 2010 and 2008 model is that the 2010 model has cor- 
rected the lifetimes of low-mass (old) TP-AGB stars. These his- 
tograms show the effect of this correction, there are very few low- 
mass (old) luminous TP-AGB stars in the 2010 models. With this 
correction, about 65% of the 2010 models are in rough agreement 
with the actual TP-AGB star counts in the observed data. 

1.6 fiui (i.e. brighter than the TRGB). To test if there 
are large numbers of NIR luminous AGB stars missing 
from our samples, we return to the NIR only CMDs 
{FllOW - i^ieOW^). We find that there are typically 
several red TP-AGBs with luminosities brighter than the 
TRGB that were missed in the optical-IR CMD search. 
However the total numbers of missing stars are typically 
fewer than the Poisson uncertainties of the original count. 
In addition, because these are among the most dust ob- 
scured, they tend to not be among the most luminous 
TP-AGBs at 1.6 ^m. Therefore missing AGB stars in 
the data cannot account for the model/data differences 
in the NIR flux. 

What about the role of dust obscuration on the model 
itself? Real AGB stars can experience significant self- 
induced dust obscur ation and drop ou t of both opti- 
cal and NIR CMDs (jBover et al.ll2009f ). In fact, some 
model prescriptions predict that for progenitor ages of 
100 - 200 Myrs, TP-AGB stars can spend nearly half of 
their lives as highly dust obscured objects ()Marigo et al.l 
l2008f ). However, these dust prescriptions have not been 
included in the SPS codes we are using to model the 
CMD. Therefore we expect that some of the model TP- 



AGB stars would be much fainter at NIR wavelengths if 
the proper dust prescriptions were included. 

One way to explore this issue further is to track the 
carbon-to-oxygen (C/0) ratio in the model AGB pop- 
ulations. Over the course of an AGB star's lifetime, 
carbon is dredged up from the interior to the surface, 
changing the overall C/0 ratio. When this ratio ex- 
ceeds unity the star is termed a carbon star. The ad- 
ditional carbon makes it much easier to form dust in 
the stellar atmosphere. Therefore the reddest TP-AGB 
stars tend to be asso ciated wi th carbon-rich popula- 
tions (e.g.l Nikolacv fc Weinberg! l 2000H Cioni et"alll2006l: 
iBonanos et al...2009; .Bover et al.ll2011[ ). 

We use the F AKE routine (which is part of the 
MATCH package iDolphinI 120021 ) to track the C/0 ra- 
tio, metallicity, age, and mass of the artificial TP-AGB 
stars. What we find is that the most luminous model 
TP-AGB stars, which are typically quite blue, are metal- 
poor and carbon-rich. Compared to the metal- rich AGB, 
metal-poor stars typically need less carbon dredge up to 
reach a C/0 ratio greater than unity. However, because 
the SPS models we are using do not include prescriptions 
for circumstellar dust, the colors of these artificial carbon 
stars appear to be driven by metallicity, rather than their 
C/0 ratio. Being metal poor, these TP-AGB stars are 
blue, even though they should have significant amounts 
of circumstellar dust. In contrast, the reddest TP-AGB 
model stars, especially in galaxies DD082, IC2574-SGS, 
and M81, are oxygen-rich stars with high (solar) metal- 
licity. High metallicity means that more carbon needs to 
be dredged up to to become a carbon star, so these very 
red TP-AGB stars are actually modeled as oxygen-rich. 
Again, their red color is driven by their metallicity rather 
than their C/0 ratio. 

These results are counter to the observationally-driven 
expectations that the reddest stars will be carbon-stars 
and the bluer ones will be oxygen-rich. A more complete 
treatment of self-obscuration by dust will be explored 
further in the next paper in this series, and may be key 
for accurate modeling of the TP-AGB, even in the NIR. 

4.1.2. Are the Measured SFHs of the Sample Accurate? 

If the input SFHs are incorrect, the SPS models are 
unlikely to match the data. We now test if the observed 
discrepancies between the model and real AGB stars can 
be explained by uncertainties in the SFHs of our galax- 
ies. To do so, we systematically lower the SFRs for inter- 
mediate aged populations (< 2 Gyrs) to bring them in 
line with the predicted numbers and fluxes of TP-AGB 
stars and then recalculate the model CMDs. We find 
that by systematically lowering the SFRs by roughly the 
uncertainties in a given age metallicity bin, we can re- 
duce the model predicted AGB contributions by roughly 
a factor of 2. While this solution appears to fix much 
of the data/model discrepancies on the AGB, it actually 
creates a larger problem in another region of the CMD, 
namely the main sequence turn-off (MSTO) for interme- 
diate aged populations. The systematically lower SFRs 
now under-predict the numbers of stars on the MSTO by 
a factor of two. The number of TP-AGB stars and inter- 
mediate aged MS stars are therefore in tension. However, 
the MSTO is much better understood than the AGB, and 
is better populated making it more robust to stochastic 
fluctuations. It is therefore more likely that the original 
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Figure 10. CMDs of three sample galaxies, selected to have had significant recent star formation, compared to SPS models (based on 
the 2010 Padova isochrones), with the models divided into different time bins. Colored regions are same as Figure [2| only now we add an 
additional luminosity cut on the TP-AGB bin to compare the numbers of the most luminous TP-AGB stars with the model predictions. 
These high luminosity TP-AGB stars are primarily at young ages (< 1 Gyr). Typically the models contain more high luminosity TP-AGB 
stars than the data. 



SFHs were reasonable. 

In addition to forcing a SFH with a systematically 
lower SFR at young ages, we also have run MC sim- 
ulations that span the full range of input SFHs deemed 
acceptable by the SPS modeling routine CalcSFH. As de- 
scribed previously, we use the MC simulations to derive 
the true uncertainties on the numbers and fluxes from 
the TP-AGB (see Figure SJ. While the uncertainties on 
any one galaxy are large (e.g. 20-70%), they are typically 
not large enough to account for the offset in flux from the 
data, especially for galaxies with significant recent star 
formation. In all cases, the best-fit SFH systematically 
over-predicts the TP-AGB flux. The weighted average of 
the model over-prediction in flux is 230% with a standard 
deviation of 80% for the ensemble of galaxies. 

4.2. RHeB Model and Data Discrepancies 

Our analysis in §3.4 shows that RHeB stars are also a 
significant contributor to the NIR flux in galaxies with 
ongoing star formation. However, the best-fit SPS mod- 
els do not capture the properties of RHeB populations 



in multiple ways. First, the SPS models tend to under- 
predict the numbers of RHeB stars (see Table 4) . Second, 
they tend to under-predict the luminosities of individual 
RHeB stars. This result can be easily seen in Figure fTOl 
where the most luminous RHeB stars in IC2574, for in- 
stance, are almost a magnitude brighter than the most 
luminous model RHeB stars. Third, the model RHeBs 
appear to be redder than the data in these NIR filters. 
This result can also be seen in Figure [101 While the ob- 
served RHeBs in IC2574 are clearly separated from the 
red edge of the RHeB selection region, the model RHeB 
stars tend to hug the red edge of the of this region, with 
a color offset of ~ 0.2 mags. Finally the RHeB branches 
in the observational data tend to be tight sequences and 
thus are easily captured by the narrow boxes in CMD 
space. The model RHeB sequences, on the other hand, 
tend to be spread out in color, forming much looser se- 
quences, which may contribute to stars falling out of the 
same boxes in the CMD. 

Figure [10] also reveals several additional issues, in- 
cluding the problematic effects of photometric uncertain- 
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Figure 11. The fraction of NIR ligtit arising from IR, luminous stars as a function of redshift as predicted by the 2008 (and 2010) Padova 
isochrones. We generate model populations at high redshift with two different star formation histories. Star formation begins at ^ = 6 
for both models, however one undergoes an extended burst that peaks at z = 3 with a SFR of 30 M_0 yr~l (left), and the other has a 
constant SFR of 5 M_0 yr~l (right). Top panels show the SFHs, while bottom panels show the predicted flux contributions from, RHeBs 
(blue), TP-AGB stars (red) and additional MS stars brighter than the TRGB (stars). While RHeBs dominate at early times (producing 
40% of the light ai z = 5.5) , TP-AGB stars dominate by z = 4, reaching a peak of almost 70% of the Fl&OW flux in the extended burst 
model and over 50% for the constant SFR. These models are representative of the types of stellar population synthesis currently used when 
fitting high-2 CMDs. However, they may need to be corrected to account for the observational results presented in this paper. 



ties when selecting subregions of the CMD. In galax- 
ies IC2754-SGS and UGC4305-1, a significant number of 
stars in the 0.3 — 1 Gyr time bin scatter into the RHcB 
selection region, even though they are likely to be AGB 
stars. This effect also drives the large model/data num- 
ber ratio for RHeBs in galaxy NGC 2976, where there is 
only one real RHeB star but a larger number of predicted 
stars. However, removing these spurious stars from the 
RHeB classification would only tend to strengthen our 
conclusion that the models tend to under-predict the 
RHeB flux in these galaxy. Of course, there may be a 
similar contribution of TP-AGB to the RHeB bin in the 
actual data, in which case no correction would be neces- 
sary. 

An even more important issue with the RHeB mod- 
els is that some model stars may be sufficiently red to 
leave the RHeB selection region altogether and instead 
fall into the TP-AGB selection region. For instance, the 
combined numbers of model RHeB and TP-AGB stars 
in the earliest time bin of Figure [10] is very similar to 



the total number of RHeBs in the data, in all three cases 
shown here. If these young red stars actually belong on 
the RHeB sequence, but are being added to the TP-AGB 
sequence, then the effect will be to over-predict the TP- 
AGB flux, and under-predict the RHeB flux, a trend that 
exists in the data. 

As with TP-AGB results, these data/model differences 
could be the result of a wide range of issues: (1) there 
could be a population of foreground stars contributing 
to the observed data, that do not exist in the model; (2) 
there could be problems with measuring accurate SFHs 
or correctly populating CMDs over the short timescales 
of RHeBs; or (3) there could be lingering issues with 
the stellar evolution codes, for example, the lifetimes of 
the most massive and luminous RHeBs could be under- 
predicted. The first issue was examined in detail in Sec- 
tion 3.6, where we found that foreground contamination 
could likely bias the observations when only a few RHeBs 
are observed. However, foregound contamination was 
found to be small compared to the observations of well- 
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Figure 12. Same as Figure 111! only now wc include corrections to the model predictions for TP-AGB and RHeB fluxes. RHeB flux 
contributions are increased by a factor of two, while the TP-AGB contributions are reduced by a similar factor compared to the 2 010 
Padova models. The over-all contribution of massive stars to the F160W flux is only reduced slightly from a peak of 80% in Figure [TT] 
to a peak of 70% here. However the makeup of the population has changed dramatically. Now the RHeB contribution dominates or is 
comparable to the TP-AGB contribution for much of the observed time-period, except for late times in the extended burst model. This 
has signiflcant implications for SEDs of galaxies at high redshift. While TP-AGB stars preferentially increase the NIR-to-optic al flux ratio 
of a p opulation, RHeBs are significantly bluer. In addition, RHeB are often accompanied by similar numbers of blue HeBs (M cQuinn et al.l 
I2OIII) . Therefore RHeBs and their blue counterparts are likely to increase both the optical and NIR luminosities of galaxies, and lower 
their M/L ratios across the full SED. Also shown are the NIR M/L ratios of these model galaxies (bottom). The M/L ratio varies rapidly 
by over a factor of 10 between z = 5 and 2 for the extended burst model (and by roughly a factor of ~ 7 for a constant SFR). To emphasize 
how these populations can rapidly change the M/L ratio of a population we now also show a multi- burst model (right). In the multi-burst 
model the M/L ratio oscillates by a factor of two as the bursts turn on and off. The amplitude of this oscillation would be even larger for 
larger burst strengths, or longer quiescent periods between bursts. Note: the widths of the spikes in the sand diagram for the multi-burst 
model (right, middle), are artificially wide because of the course redshift sampling. 



populated RHeB sequences. We examine potential issues 
with modeling young SFHs below but leave any updates 
to the stellar evolution code for the next paper in this 
sequence Rosenfield et al. (in preparation). 

4.2.1. Are the SFHs and SPS Codes Accurate for the 
Youngest Populations? 

Two of the most obvious issues with the RHeB 
branches are that the best-fit SPS models are too diffuse 
in color space and do not contain stars bright enough 
to match the data. Both of these differences could po- 
tentially be caused by inaccurate SFHs at these young 
ages. For instance, our best-fit SFHs have a relatively 
broad metallicity spread (typically 0.25 dex) even at 
young ages. This metallicity spread will cause the mod- 
eled RHeB branch to be broader than if the metallicity 
spread were narrower. We estimate that dropping the 
metallicity spread to dex could increase the numbers 
of synthetic RHeB stars by as much as 25%. While not 
enough to account for a factor of 2 discrepancy between 



the models and the data, this change would tend to im- 
prove the match between the models and data. To test 
this possibility more rigorously, we recalculate the opti- 
cally derived best-fit SFHs with CalcSFH only now we 
force the narrowest possible spread on the input metal- 
licities at young ages (using the -zinc fiag, giving an ef- 
fective spread in metallicity of 0.1 dex). The SPS models 
based on these new SFHs still under-predict the RHeB 
flux by roughly the same factor as before. 

Another issue is that the SPS models do not contain 
enough very luminous RHeB stars. The stellar evolution 
models can make stars that are as bright as the data, 
but only at the youngest ages, e.g. < 50 Myrs. It may 
therefore be that the best-fit SFHs are systematically 
under-predicting the SFRs of our galaxies at the youngest 
ages. This bias would tend to decrease the numbers of 
the most luminous RHeBs and significantly decrease the 
predicted RHeB fiux contribution. 

To test if the SFHs are wrong at very young ages, 
we rerun the SFHs only now increasing the SFRs at 
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the youngest ages. Doubling the SFRs at these ages is 
enough to roughly match the actual RHeB flux observed 
in the real galaxies. Unfortunately, doing so increases the 
numbers of very young main sequence and main sequence 
turn-off stars as well. Compared to the data, these new 
models have roughly five times the numbers of stars in 
these additional regions of the CMD. Thus it is unlikely 
that SFH problems alone can account for the model/data 
discrepancies. 

4.3. Implications for Observations at High Redshift 

Our results demonstrate that short lived TP-AGB and 
RHeB stars can contribute significant fractions of the 
NIR light of local galaxies, while contributing negligible 
amounts of stellar mass. We expect that the fiux con- 
tributions of these stars will be much greater at high- 
redshift where star formation rates are high, and the 
RGB is less developed. In this section we attempt to 
put some constraint on the contribution of TP-AGB and 
RHeB stars to the rest-frame NIR fluxes of galaxies at 
high redshift. We again focus on the rest-frame NIR 
fluxes where these stars are most luminous and where 
rest-frame 1.6 /im fluxes of high-z galaxies are well con- 
strained from deep Spitzer imaging in the 3 — 8 /xm 
observed-frame bands. This wavelength regime is also 
where the future James Webb Space Telescope will be 
operating. 

First we show the predictions from the 2008/2010 
Padova isochrones, as these, and similar isochrones, are 
the foundation for interpreting colors and magnitudes at 
high-z. Then, by reducing the contribution of TP-AGB 
flux and increasing the contribution of RHeB flux, we 
produce predictions that could account for the discrep- 
ancies between the SPS models and our observational 
datasets. Finally, we discuss these predictions in the con- 
text of SED fitting at high redshift. 

We caution that our observational constraints at low 
redshift are primarily tracing low metallicity systems, 
typically 1/10 - 1/5 solar. Massive high redshift sys- 
tems such as sub- mm galaxies and BzK g alaxies are likely 
to be at higher mean met allicity (e.g.. iSwinbank et al.l 
120041 : lOnodera et al.l l2010( ) , despite the over all decline 
in gas phase metallicity wi th redshift (e.g., lErb et ahl 
120061 : [Moustakas et al.l [20Tll) . However, our results are 
directly applicable to lower mass systems at high red- 
shift such as z = 3 Lyman break galaxies ("Er b et aD 
|2006; Ma nnucci et a l. 2009; So mmariva et al. ,20lll) and 
gamma ray burst host systems ( Laskar et al.ll20Tll ) which 
have metallicities 1/10 ~ 1/2 solar. Our results will also 
be applicable to the even lower mass systems that will 
be observed with the James Webb Space Telescope and 
the Thirty Meter Telescope. 

We use CalcSFH and NoisyCMD to generate model 
stellar populations at high redshift with two different 
star formation histories, an extended burst and a con- 
stant star formation history. We assume that star forma- 
tion begins at z = 6 for both models, and we track the 
model galaxies through z = 1, a time- frame of roughly 5 
Gyrs. The extended burst reaches a maximum star for- 
mation rate of 30 Mq yr~^ at z — 3, and then declines 
to a SFR = Mq yr~-^ by z = 1. The constant star 
formation model has SFR — 5 Mq yr~^ over the full 
time period. These star formation rates are 2 — 4 orders 
of magnitude higher than those measured in the dwarf 



galaxies in our study, but they are representative of the 
SFRs found in the high-z galaxies studied with current 
instrumentation. In each case we assume that metallicity 
increases from [M/H] = —0.7 to by z = 1. We create 5 
burst and 5 constant SFR models at each redshift, and 
average the results for each set. 

FigurelTTjshows the predicted fraction of the NIR lumi- 
nosity arising from NIR luminous stars for high-z galax- 
ies modeled with the Padova 2008 isochrones (We use 
the Padova 2008 models here as they are most similar to 
the AGB-heavy models that have been used to date in 
interpreting high-z data. However, the 2010 isochrones 
will give roughly equivalent results because there has not 
been enough time to create significant numbers of low- 
mass TP-AGB stars). When star formation commences, 
RHeB stars rapidly become the dominant IR source in 
our model galaxies, providing ~ 40% of the F\QQW light 
at z = 5.5. At that time, TP-AGB stars and mas- 
sive main-sequence stars provide another small fraction 
(~ 10%). The remainder, ~ 50%, comes from less mas- 
sive main sequence stars. By z = 3.5, the contribution 
from RHeBs declines due to increasing numbers of lumi- 
nous intermediate-aged TP-AGB stars. At this redshift, 
TP-AGB stars are predicted to be the dominant NIR 
source, producing 50% or more of the total. 

Figure [TT] gives an overview of the current predictions 
for massive-star fiux contributions to high-z galaxies. 
These fractions are baked into current high-z SED fit- 
ting and mass estimates. However, our data suggest that 
there may be problems with these predictions. In par- 
ticular, the RHeBs contribution may be under-predicted 
by as much as a factor of two, and the TP-AGB fraction 
may being over-predicted by a similar factor. Taking 
these approximate corrections into account, a different 
picture emerges for the luminosities of high-z galaxies, 
as shown in Figure [121 Now the RHeB phase dominates 
or is comparable to the TP-AGB phase, except for late 
times in the extended burst model. 

Coincidentally, there is only a small change in the over- 
all flux of NIR light from massive stars between Figures 
[m and [H The dechne in TP-AGB flux is effectively 
offset by the increase in RHeB flux, at least for these rel- 
atively short SFHs. The NIR M/L ratios are therefore 
not radically different between the 2008 models and the 
approximately corrected version. However, the different 
make-up of the population has significant ramifications 
for SED fitting. TP-AGB stars will tend to increase the 
NIR-to-optical flux ratio, because TP-AGB stars are very 
red (more luminous in the NIR than the optical) . In con- 
trast, RHeB stars are significantly less red than TP-AGB 
stars. In addition, the existence of RHeB stars usually 
implies the existence of Blue Helium Burning (BHeB) 
stars, which are very luminous a t bluer wavelengths (e.g. 
iDohm-Palmer fc Skillmanl 120021 ). Massive core helium 
burning stars tend to evolve through blue and red loops 
before ending their fusion lifetimes, and the numbers of 
BHcB are typically co mparable to the numbe rs of RHcBs 
at a given luminosity ([McOuinn et al.l 120111 ). Thus to- 
gether the red and blue core helium burning populations 
will tend to simultaneously increase both the optical and 
the NIR luminosity rather than affecting the NIR alone 
as expected for TP-AGB stars. 

These results may help explain several puzzling out- 
comes for SED fitting at intermediate redshift. For in- 
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stance iKriek et al.1 (|2010D found that SEDs of z - 0.7 
post-starburst galaxies cannot be simult aneously fit in 
the optical and NIR with lMarastonI (I2005D models , which 
have a large TP-AGB component. The iMarastonl (|2005l ) 
models tend to favor larger NIR fluxes for a given op- 
tical flux, compared with actual observations. However, 
both the blue and r e d side s of these SEDs can be fit with 
iBruzual k, Charlotl (|2003D models, which have a s maller 
TP-AGB contribution. Likewise, iMuzzin et all H2009f ) 
found that SEDs of massive compact galaxies at z ~ 2 
are better fit with Bruzual & Chariot (2003) models than 
updated versions that include a larger TP-AGB compo- 
nent. In both cases the AGB-lite models tend to predict 
more stellar mass than the AGB-heavy models. 

Our results also suggest that AGB-lite models may be 
a better match to the data from lo cal galaxies . This 
alone may be enough to explain the IKriek et al.l ()2010l ) 
and lMuzzin et al. (2009) results. However, their conclu- 
sion, that these post-starburst galaxies contain higher 
stellar mass than the AGB heavy models would predict, 
may be incorrect. If their sample galaxies contain some 
stellar populations younger than ~ 300 Myrs, then red 
and blue HeB stars may also be contributing flux. By 
increasing the fiux at both the optical and NIR, these he- 
lium burning stars necessitate a lower M/L ratio across 
all wavelengths. The current AGB-lite models do not 
necessarily adequately capture the HeBs, and our data 
suggest that they must be considered when there has 
been recent star formation. 

Figure [12] also shows the evolution of the NIR M/L 
ratio for our example high-z SFHs (the extended burst 
model and the constant SFR model) after applying ap- 
proximate corrections that increase the RHeB fiux and 
decrease the TP-AGB fiux. For both SFHs, the M/L ra- 
tio changes rapidly, and by large factors — over a factor 
of 10 between z = 5 to 2 for the extended burst and a fac- 
tor of 7 for the constant SFR. Capturing this rapid evo- 
lution is essential for estimating accurate stellar masses 
and proper SED fitting at high-z. For example. Figure 
[T2l shows a multiple-burst SFH, where the M/L ratio os- 
cillates by a factor of roughly 2 as the bursts turn on and 
off. These oscillations amplitudes will be even stronger if 
the burst strengths are larger or the quiescent intervals 
are longer. 

Again, we note that Figure [12] is intended to demon- 
strate the consequences of the largest reasonable changes 
to the stellar evolution codes, based on our observational 
data. If some of the model/data differences described in 
the previous sections arise from poor SFH measurements, 
or issues with the SPS modeling, then the changes we 
are describing here will be smaller than we are report- 
ing. A final assessment of the model/data comparison, 
and model updates will be provided in future papers in 
this series. 

5. CONCLUSIONS 

We report on the contribution of evolved stars to the 
NIR fluxes of 23 nearby galaxies observed with HST 
WFC3 (NIR) and ACS (optical). The HST color- 
magnitude diagrams separate different phases of stellar 
evolution, including the red giant branch, the asymptotic 
giant branch and luminous red helium burning stars. The 
CMDs also provide markers of the star formation histo- 
ries of each galaxy. We use the optically derived SFHs to 



model the NIR CMD to the main sequence K dwarfs. We 
use these models to produce an estimate of the total 1.6 
/im flux for each galaxy by combining the observed NIR 
fluxes for luminous stars with the model fluxes of fainter 
stars. We then calculate the contribution of TP-AGB 
and RHeB stars to this total flux. 

While the RHeB and TP-AGB sequences represent 
negligible stellar mass, they can account for as much as 
21% (or 18% after foreground correction) and 17% (re- 
spectively) of the 1.6 /im fluxes for these local galaxies. 
At higher redshift, when galaxies do not have well de- 
veloped red giant branches, RHeBs and TP-AGB stars 
are expected to produce much higher fractions of the 
NIR light. At high-z, the summation of all IR luminous 
sources (the rare stars brighter than the typical TRGB at 
that time) can be expected to be as high as 70%, most of 
which comes from TP-AGB and RHeBs. These massive 
stars can therefore produce short-lived dramatic changes 
in the NIR M/L ratio. Stellar masses based on NIR 
fluxes could therefore be fraught with large systematic 
uncertainties. Likewise, these evolved populations may 
contribute to the scatter of NIR galaxy scaling relations 
such as the NIR metallicity-luminosity and Tully-Fisher 
relations. 

We compared our observational results to predic- 
tions from optically derived, best- fi t SPS models ( based 
on the isochrones from Padova, iMarigo et al.l 120081: 
iGirar di et al. 2010). We confirmed the iGirardi et al.l 
(j2010t) resuh that SPS models based on the 2008 Padova 
isochrones over-produce the old (low-mass) TP-AGB 
populations. We also showed that SPS models based on 
the 2010 version of the Padova models have largely cor- 
rected the TP-AGB numbers problem; the SPS models of 
65% of our sample match the observed numbers within 
the measured uncertainties. While a small systematic 
in the numbers remains (Nmodei/Ndata = 1-5 ± 0.5, or 
1.4 ± 0.5 after removing outliers) a larger bias exists in 
the predicted FIQQW fluxes from TP-AGB populations, 
with Lynodei/ Ldata ^ 2.3±0.8 (or 2.2±0.8 after removing 
outliers). 

This model/data discrepancy in TP-AGB flux is pri- 
marily driven by predictions of very luminous (and 
young) TP-AGB stars that are largely absent from the 
NIR data. We show that uncertainties in the measured 
SFHs are unlikely to account for the bulk of this over- 
prediction. Part of the discrepancy is likely due to the 
fact that galaxy simulations in this paper do not take into 
account the TP-AGB pulse cycle luminosity variations 
and, more importantly, obscuration of the optical and 
NIR due to dusty circumstellar envelopes. Since both 
aspects are included in the original datasets of evolution- 
ary tracks and isochrones, their effect will be explored in 
the next paper on this series. 

This work also suggests that the typical uncertainties 
of < 50 % in the lifetimes of TP-AGB models, which seem 
to have been reached by present-day calibrations based 
on Magellanic Cloud data, may still not be good enough 
for the accurate modeling of distance galaxies and the 
derivation of their stellar masses, because of the possible 
large contribution of a few bright TP-AGB stars to the 
integrated light. If our next goal is to reduce typical 
errors in the lifetimes and fluxes of TP-AGB populations 
by a factor of at least two, for a wide enough range of 
stellar masses and metallicities, it is clear that such a 
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goal cannot be reached using Magellanic Cloud data only. 
It requires large samples of nearby galaxies with well- 
constrained SFHs as those presented in, e.g., Gullieuszik 
et al. (2008b), Melbourne et al. (2010), Girardi et al. 
(2010), and in this paper. 

We also examine the predictions of the best-fit SPS 
models for RHeB stars. In general they predict fewer, 
fainter, and redder RHeB stars than are observed in our 
galaxies. As a result the models typically under-predict 
the flux contribution of RHeBs by a factor of 2.0 ± 0.6 at 
F160W, for galaxies with significant RHeB populations. 
As with the TP-AGB stars, issues with the measured 
SFHs are unlikely to account for all of the data/model 
discrepancies on the RHeB branch. However, lingering 
issues with how the SPS models populate the CMDs or 
with the stellar evolution isochrones themselves could re- 
solve these differences. 

If we take our mo del/ data differences at face value, 
we can make approximate predictions for how the fu- 
ture models might behave when used to study more dis- 
tant galaxies. For instance current AGB-heavy models 
predict higher rest-frame NIR-to-optical flux ratios than 
our data would support, a discrepancy ha s been seen in 
several high-z SE P fitting programs (e.g. iMuzzin et al.l 
I2009t iKriek et al.|[20ia) . These results will be explored 
further in Rosenfield et al. (in preparation), and, if 
needed, any updates to the Padova isochrones will be 
included at that time. 
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